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Abstract 

Quantum Field Theory on Curved Spacetime (QFT on CS) is a well established theoretical framework which intuitively 
should be a an extremely effective description of the quantum nature of matter when propagating on a given background 
spacetime. If one wants to take care of backreaction effects, then a theory of quantum gravity is needed. It is now widely 
believed that such a theory should be formulated in a non - perturbative and therefore background independent fashion. 
Hence, it is a priori a puzzle how a background dependent QFT on CS should emerge as a semiclassical limit out of a 
background independent quantum gravity theory. 

In this article we point out that the Born - Oppenheimer decomposition (BOD) of the Hilbert space is ideally suited in 
order to establishing such a link, provided that the Hilbert space representation of the gravitational field algebra satisfies 
an important condition. If the condition is satisfied, then the framework of QFT on CS can be, in a certain sense, 
embedded into a theory of quantum gravity. 

The unique representation of the holonomy - flux algebra underlying Loop Quantum Gravity (LQG) violates that 
condition. While it is conceivable that the condition on the representation can be relaxed, for convenience in this article 
we consider a new classical gravitational field algebra and a Hilbert space representation of its restriction to an algebraic 
graph for which the condition is satisfied. An important question that remains and for which we have only partial answers 
is how to construct eigenstates of the full gravity - matter Hamiltonian whose BOD is confined to a small neighbourhood 
of a physically interesting vacuum spacetime. 
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1 Introduction 



Quantum Field Theory (QFT) on Curved Spacetime (CS) is nowadays a well established discipline of mathematical 
physics (see [T] for a recent review and [5] for a functorial formulation) . It is designed to describe the regime of physics 
in which backreaction effects of matter on geometry and the quantum fluctuations of geometry itself can be safely 
neglected. In order to go beyond that, a theory of quantum gravity is needed. It is nowadays widely believed that 
perturbative quantum gravity is non renormalisable and therefore a non perturbative formulation has to be found. Of 
course, such a theory of quantum gravity must contain QFT on CS in its semiclassical limit. The puzzle that poses 
itself is, how a background independent theory can describe a background dependent theory such as QFT on CS all. 
This question is the central topic of the present paper. 

We still do not have a working theory of quantum gravity, however, there are several proposals which are still 
under construction. One of them, Loop Quantum Gravity (LQG) PUS] indeed is a background independent approach 
to the quantisation of general relativity and in the last 20 years a lot of effort has been put in the rigorous definition of 
the mathematical framework of this theory. Background independence is a feature that distinguishes LQG from most 
other quantum field theories constructed so far, as for example the standard model of elementary particle physics. 
When examining ordinary quantum field theory with care one realises that much of the structure we are used to from 
the particle physics point of view, including for instance the construction of a Fock space as the fundamental Hilbert 
space of the theory, the concept of a particle as a certain mode-excitation of the field, the notion of causality and 
the notion of a vacuum state, are strongly linked to the presence of a given background spacetime. For a quantum 
theory of general relativity there is no such background spacetime available because we are quantising the dynamics of 
spacetime itself, and thus it should not be too surprising that the fundamental structures appearing in a quantisation 
of general relativity differ a lot from the ones used in ordinary quantum field theory. 

The natural Hilbert space to describe the quantum dynamics of LQG is not a Fock space as one might expect - 
in fact the Hilbert Hail space used in LQG, named the Ashtekar-Isham-Lewandowski Hilbert space (5j E] , is the 
space of square-integrable functions over a set of appropriately generalised SU (2)-connections on a 3 dimensional 
manifold and does not bear a lot of similarity with the Hilbert spaces known from ordinary QFT. However, when 
trying to find irreducible representations of the kinematical algebra of GR, also called the holonomy - flux - algebra 
in the LQG context, one is naturally led to this space. Analogous to the well known Stone-von Neumann theorem 
in quantum mechanics, there exists even an uniqueness result [8] , that states that under certain assumptions ( 
spatial diffcomorphism invariance being the main one) the representation used in LQG associated with 7Yal is unique. 
The structure of this Hilbert space is quite different from the ones known from ordinary quantum field theory on a 
given background: It is manifestly background independent and its elements can be written in terms of so called spin 
network functions, functions living on certain embedded graphs which are labeled by irreducible representations of 
SU(2). Furthermore, the flux operators, corresponding to the momentum operators in the canonical picture, encode 
geometrical information about the quantum state turn out to fulfil a non-commutative algebra. 

The relation between LQG and classical general relativity has been studied in detail [HI HH1 [HI HH HHI and it turns 
out that LQG, or rather its recently proposed modification coined algebraic quantum gravity, see [H| H31 Hj ) , has the 
right semiclassical limit at least as far as the infinitesimal dynamics is concerned. 

However, if LQG is a fundamental quantum theory for gravity and matter, then it should be possible to rediscover 
ordinary quantum field theory on a fixed background within a sector of the fully background independent framework 
of LQG. There has been some work into this direction |16[ I17j. but some issues remained unresolved: First, one did 
quantise at the level of the kinematical Hilbert space. Therefore all operators in question were not gauge invariant, the 
coherent states that one was using were not gauge invariant and thus one might rightfully ask whether the constructions 
displayed there would survive upon passage to the physical Hilbert space. Next, the matter Hamiltonians used there 
were just the matter contributions to the Hamiltonian constraint and one did not have a convincing argument why 
those should be used as Hamiltonians driving the physical time evolution of the system. Finally, backreaction effects 
were not taken into account. 

The aim of this article is to try to improve the understanding of how a quantum field theory of matter on a given 
background can emerge out of a theory of quantum gravity. In order to deal with the afore mentioned constraints we 
carry out a reduced phase space quantisation of general relativity plus matter degrees of freedom using a dynamical 
reference frame given by some dust fields using methods from [18, 19, 20 . Thus we arrive at a true (i.e. non-vanishing) 
physical Hamiltonian that generates time evolution in the reference frame described by these dust fields. Since in 
contrast to the total Hamiltonian constraint the physical Hamiltonian is non - vanishing, in general, the physical 
Hamiltonian operator will involve more complicated terms than just the total Hamiltonian constraint. In particular 
it involves a square root of polynomials of the contributions to the constraints associated with the non - dust degrees 
of freedom. To deal with this square root is a difficult technical subject of its own which is not the main concern 
of this paper. Thus, we consider a regime of the theory for which the Taylor expansion of the square root can be 
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restricted to the dominant term which is just the integrated non - dust contributions to the Hamiltonian constraint. 
Therefore we get an a posteriori justification for why one can build a matter QFT just using the matter parts of the 
Hamiltonian constraints as dynamical operators which otherwise would seem rather ad hoc. Of course, there is also 
the gravitational contribution which we will take care of in this paper. That this is possible relies on the judicious 
choice of the Brown - Kuchaf dust fields as rods and clocks. 

When examining LQG and ordinary QFT from a mathematical point of view it turns out that there are two main 
issues which need to be clarified when the latter should emerge out of the former in some appropriate limit: First, 
although LQG is a true continuum theory, its basic building blocks are defined on graph like structures. On the other 
hand, the Fock space of ordinary QFT uses the continuum properties of the underlying spacetime and one has to 
understand how one can arrive at such a space starting with a discrete theory at least in some limit. 
Second, the flux operators in LQG, which are the basic building blocks if we were to construct a 'metric operator' as 
a quantum analog of the classical metric form a non-abelian operator algebra in LQG. This non-commutativity origi- 
nates in the fact that the fluxes are obtained by smearing the sensitised triads along hypersurfaces of codimension one 
(and thus, singular ones from the three dimensional point of view) |21j which causes the associated Hamiltonian vector 
fields to be non commutative. Hence in any representation of this algebra on a Hilbert space, the non commutativity 
prohibits that the three geometry operators which couple to the matter fields can be simultaneously diagonalised. 
This implies an obstacle to directly applying the framework of QFT on CS when choosing a representation for the 
matter fields because there one relies on a commutative background spacetime. While much work on QFT on non 
commutative spacetimes has been conducted (see e.g. [HJ H2 HH HH] and references therein), there one assumes a 
non commutative structure in the coordinates rather than the metric and therefore we are not in the position to apply 
that framework. Thus, if one wishes to import QFT on CS techniques into LQG, one first has to understand how to 
"Abclianise" the metric operators in appropriate regime. 

In this article we will focus on the first point and leave the latter open for future research: Namely we will con- 
struct a quantum theory for gravity plus matter that is close to LQG in the sense that it is defined on a fundamentally 
discrete entities (algebraic graphs), but the geometrical operators are defined as ordinary (commuting) multiplication 
operators. This is possible by choosing an algebra different from the holonomy - flux algebra which shares some of 
its features (namely the discreteness of its fundamental building blocks) but deviates from it in other essential points 
(mainly the non-commutativity of geometrical operators is avoided). We choose a representation of this algebra which 
is very close to a (background independent) Fock representation^ and therefore the discreteness of the spectrum of 
geometrical operators of LQG is not true in this representation^ Using this algebra we examine how a Fock space for 
the matter degrees of freedom on a fixed background can emerge out of a full quantum theory for gravity plus matter. 

The idea then for how to import QFT on CS techniques into quantum gravity is to employ a Born-Oppcnhcimer 
Decomposition (BOD) of the Hilbert space: In the chosen representation for the gravitational field, the wave functions 
depend on the co-triad configuration e. For each co - triad configuration we consider a matter Hilbert space H. e labelled 
by the co - triad. That matter Hilbert space is naturally chosen as a Fock space defined by the natural creation and 
annihilation operators on the formally ultrastaticH spacetime defined by the co - triad, whenever it is non - degenerate. 
When the co - triad is degenerate, the QFT on CS construction does not work and in fact the Hamiltonian operator 
is ill - defined in this case. Its domain of definition therefore must be chosen as to exclude this case. The fibre Hilbert 
spaces H e are then glued together in a direct integral with respect to the spectral measure defined by the projection 
valued measures of the commutative algebra of co - triad operators. The cotriad operator and the momentum operators 
respectively act by multiplication and (suitable) derivation with respect to e in the fibres respectively. Also the matter 
field operators, being expressible as linear combinations (with e dependent "coefficients") of creation and annihilation 
operators, preserve these. This works when the Hilbert space for the total system is separable. Conversely, choosing 
a spectral measure and an action of the momenta conjugate to the triads with respect to which they are self - adjoint 
defines a separable Hilbert space representation. Notice that in this construction the question of unitary equivalence of 
the fibres Tt e never arises (they will typically be unitarily inequivalent). We will generalise this construction and start 
with a non - separable Hilbert space based on von Neumann's Infinite Tensor Product. Now the separable Hilbert 
spaces Ti. e (and uncountably infinitely many more) arise as strong equivalence class separable sectors in the matter 

1 Essentially, we smear the continuum co-triad along paths and the conjugate momentum along surfaces, similar as connections and 
densitised triads in LQG. Notice that the smeared co-triads in contrast to the holonomy of a connection do not transform locally under 
gauge transformations which makes it difficult to construct gauge invariant objects from them, in contrast to the Wilson loops. However, in 
this paper we consider a top to bottom approach and define an abstract algebra on an algebraic graph which approximates the continuum 
algebra and which transforms locally under gauge transformations. 

2 There is no contradiction to the uniqueness theorem because we use a different algebra. 

3 A spacetime metric is called ultrastatic if it is static (i.e. stationary and the shift vanishes) and the lapse equals unity. Thus the non 
- trivial components of the four - metric are encoded in the three - metric which is not explicitly time dependent. There is no explicitly 
time dependence here because our Hamiltonian is not explicitly time dependent by virtue of the choice of the material reference frame. 
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ITP labelled by the corresponding vacua fi e . They are automatically orthogonal in the ITP inner product when the 
fi e lie in different strong equivalence classes. The advantage of the ITP over the separable Hilbert space construction 
is that we do not need to choose a spectral measure (and the TL e for degenerate TL e which could be defined e.g. by 
an e independent Fock representation) and that it is easier to find a self - adjoint action of the operator conjugate 
to the triad. Moreover, the ITP contains vectors of the form f2 g (g> fi m where £l g ,tt m lie in the gravitational and 
matter ITP respectively. This not possible in the separable construction (at least not for non compact spacetimes for 
which an infinite algebraic graph is essential) because this would mean that the same il m lies in all the TL e which is a 
contradiction when the 7i e are not all unitarily equivalent. 

The use of Born-Oppenhcimcr Decomposition in the context of quantum gravity has a long tradition, sec for 
example [27] and references therein for a historical review of the Born-Oppenhcimer method in quantum gcomctrody- 
namics or |28| for a recent application in the context of spinfoam models. However, to the best of our knowledge, the 
connection between the BOD and QFT on CS has not been much exploited yet. As already mentioned, one would like 
to use the BOD also for the representation underlying LQG, however, at the moment there seems to be no promising 
idea for how to do this when those slow degrees of freedom that couple to the fast degrees of freedom cannot be 
simultaneously diagonaliscd. 

The structure of this article is as follows: 

In section [5] we very briefly review how one can construct a true physical Hamiltonian for general relativity plus 
arbitrary matter at the classical level using a preferred reference frame given by four dust fields. For details concern- 
ing this construction see [29l HH [19]. Then we define a regime for which it is justified to approximate the physical 
Hamiltonian by just the integrated contributions of the non - dust degrees of freedom to the Hamiltonian constraint. 
In the remainder of the article we will use this approximation to define dynamics in the dust frame. Next, we rewrite 
classical GR in terms of an 5'J7(2)-cotriad and its canonically conjugate momentum which is essentially a densitiscd 
extrinsic curvature. The smearing of these objects along one dimensional paths and two dimensional surfaces respec- 
tively defines a new Poisson algebra for classical GR which separates the points of the phase space and which we 
will call the (extrinsic) curvature - length algebra because the square root of the trace of the square of the cotriad 
along short paths is approximately the length along that pathQ. In section [3] we construct a representation of the 
curvature - length algebra and of a scalar field of Klein-Gordon type0 on an abstract algebraic graph in the spirit 
of [TJ] [T31 [23 120] and the essential feature which distinguishes this theory from ordinary LQG is that cotriads are 
quantised as multiplication operators which will be essential for the semiclassical analysis in the second half of this 
article. 

In section [4] we construct semiclassical states for the gravitational sector of this quantum theory. 
In section [5] we introduce the Born - Oppcnheimer Decomposition for our system and the associated approximation 
scheme for the full spectral problem, which still takes backreaction effects into account. Here the gravitational variables 
are treated as the 'heavy' ones and the matter variables as the 'light' ones. Such an approximation scheme is justified 
in regimes where a description in terms of interacting quantum fields on an almost classical background spacetime 
makes sense because for time scales on which particle interactions take place the change in geometry are truly small. 
We illustrate how this works using a simple minisuperspace model. In the appendix we list some facts about the ITP 
for the benefit of the reader. 

2 Classical theory 

2.1 Brown— Kuchar dust reduction and observables for General Relativity 

The problem of observables, that is, calculating quantities which are invariant under the gauge group of the theory is 
much more severe in general relativity compared to other gauge theories such as electrodynamics or SU(N) Yang-Mills 
theories. The reason is twofold: First, in general relativity the gauge group is given by Diff(M), the diffeomorphism 
group of the four dimensional underlying manifold. This leads to a rather complicated constraint algebra (often called 
the hypersurface deformation algebra or Dirac algebra in the literature) which is neither finite dimensional nor an 
honest Lie algebra (it has phase space dependent structure functions and not structure constants). Second, the issue 
of observables is strongly related to the background independent nature of the theory: Taking the manifold M just 

4 It is straightforward to construct a representation of this algebra in the continuum very close to the one of LQG by taking as 
configuration variables the exponentials of the cotriads (one for each direction in the Lie algebra) along paths which would substitute the 
holonomics of LQG and then one would use the Ashtekar — Lewandowski measure for C/(l) 3 . However, the construction of gauge invariant 
functions in the continuum from these Abelian holonomies works only in the limit of short paths which is a limit too singular in this 
representation because it is not strongly continuous under path deformations. This is why we proceed to the algebraic formulation here. 

5 A generalisation to gauge fields or fermions is straightforward but we will only consider a scalar field for conceptual clarity. 
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as an auxiliary structure, it is not meaningful to talk about fields evaluated at some point on this manifold. The true 
gauge dependent degrees of freedom should be defined as relations between dynamical fields. 

Such relational observables were introduced first in [30] [3T] 1321 [33J and later refined in [3H [35J ESI 133 [33] ■ In retrospect 
it turns out that there is a one to one correspondence between a choice of gauge fixing conditions with an associated 
reduced phase space and a choice of relational observables. The correspondence is established by choosing as rods and 
clocks of the relational observables precisely the gauge fixing conditions of the reduced framework, see e.g. [18] 39J. The 
explicit expression of these observables in terms of the non gauge invariant degrees of freedom is rather involved, see 
e.g. [101 SI] for a perturbative treatment. Fortunately, it is not needed for what matters is the Poisson algebra between 
the observables (i.e. true degrees of freedom) and the associated physical Hamiltonian (i.e. reduced Hamiltonian) [42] . 
The reduced phase space for vacuum general relativity has a very complicated symplcctic structure because the rods 
and clocks that one can construct necessarily lead to spatially non - local objects. However, in the presence of suitable 
matter these complications can be avoided. Here by suitable we mean matter that comes as close as possible to the 
mathematical idealisation of a test observer fluid and such that the physical Hamiltonian that drives the physical time 
evolution of the observables comes as close as possible to the standard model Hamiltonian [13] when the spacetime is 
close to flat. In [THEl such a suitable matter system was identified as the pressureless dust matter fields introduced 
by Brown and Kuchar in [29]. One considers the enlarged system Seh + Smatter + Sdust where S'eh is the usual 
Einstein-Hilbert action, S ma tter includes all possible other standard matter that one likes to couple to gravity and the 
dust action Sdust is given by 

S dust = -\f d A X^\d C t{g)\p[g^U^U v + 1} (2.1) 

M 

Here M is the four dimensional (spacetime) manifold which can topologically be identified with R x X for some three 
dimensional manifold X of arbitrary topology and X £ M are local coordinates on M. g^(X) with fx, v = 0, 1, 2, 3 
denotes a (Pseudo-) Riemannian metric on M. U <G T*M is a one form defined as U = —dT + WjdS J , with j = 1, 2, 3, 
for some scalar fields T,Wj,S^ <E C°°(M). So finally the action (|2.f [) is a functional of g^ v and eight scalar fields 
P,T,W h S\ 

The coupled system S'eh + S ma tter + Sdust seems to be very complicated at first sight, and indeed when performing a 
Hamiltonian analysis one realises that it is a second class constrained system. However, the second class constraints 
can be solved and on the reduced phase space the Dirac bracket is identical to the usual Poisson bracket. Moreover, if 
one chooses T, S J as the four clock fields then it turns out that one can explicitly construct observables associated to 
the spatial three metric q a b{x) and their conjugate momenta p ab (x) in the dynamical reference frame (t, <r a ) defined 
through the readings of the dust fields T, S- 7 . We will not give the details of this construction here and the interested 
reader should consult [TU [19] for all the details. What will be important for the purpose of this article is that we can 
construct a physical Hamiltonian -ff p h ys that generates physical equations of motion for the gauge invariant 3-metrics 
Qab(o~) and their respective momenta P ab (a). These observables are functions on the "dust manifold" S, i.e. ueE 
does not label a point in X but a point in the range of the clock fields S- 7 . The explicit form of the physical Hamiltonian 
reads 



H phys := I da^C 2 - Q ab C a Cb = J daH(a) (2.2) 

s s 

where C := C 9 {Q, P) + C m (Q, U, $) and C a := C°(Q, P) + C™{Q, U, $) are the contributions to the Hamiltonian and 
diffcomorphism constraints of general relativity and the matter sector written in terms of gauge invariant quantities 
Q a b(cr), P ab (a) and the corresponding matter observables denoted by n(<r), $(0-). Also we introduced the Hamiltonian 
density denoted by H(a). 

Note that C and C a are non vanishing because only the total Hamiltonian and diffcomorphism constraint including 
the contribution of the dust vanishes and the energy momentum density of the dust, which are constants of the 
physical motion, must be non vanishing anywhere in order that it provides a bona fide system of coordinates. Note 
that the matter degrees of freedom should be just understood symbolically since we have not yet decided which 
matter components wc actually want to consider. In general the algebra of observables involves a Dirac bracket 
rather than the Poisson bracket. However, for this particular choice of clock variables, Q a b and P ab as well as the 
corresponding matter observables form mutually commuting canonical pairs, that is their Poisson bracket is given by 

6 We will use capital letters to denote the gauge invariant observables in the reference frame given by the dust clocks T, S° . Thus, by 
Q ab we denote the gauge invariant three metric associated to the gauge variant 3-mctric q a i, and by P ab the gauge invariant canonical 
momenta associated to the gauge variant momenta p a t,. For the exact construction of these gauge invariant quantities see the articles cited 
above. 
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{Qabic), P cd (<r')} = 6? a S&8(a, a 1 ) and one can show that the physical equations of motion are given by 

Qab = {-Hphys, Qab\, P" b = {-^phys, P ab } (2.3) 

and likewise for the matter sector. 

The consistency of the dust picture with the usual gauge variant description of black hole and cosmological spacetimes 
has been verified in [T8l H"9l l44l . 



2.2 Approximation for the physical Hamiltonian 

The physical Hamiltonian shown in equation (|2.2j) includes a square root under the integral. This square root originates 
from solving the Hamiltonian constraint for GR coupled to dust and possibly other standard matter for the dust 
momentum P, the momentum conjugate to the dust time T. We want to discuss how the reduced framework when 
quantised is related to quantum field theory on curved spacetimes. In contrast to the standard framework we will also 
quantise the gravitational sector of the theory by quantising H p hy a which contains all (gauge -invariant) degrees of 
freedom except the dust ones. Although the physical Hamiltonian H p h ys has been quantised in its full form using Loop 
Quantum Gravity techniques and spectral theory [20] , in this paper we will consider i? p hy S in a certain approximation 
which we will now discuss. 

First of all, the derivation of -ff p h ys in [TS] shows that one has the two anholonomic constraints C 2 — Q ab C a Cb > 
and C > on the constraint surface as otherwise the momentum P conjugate to T becomes imaginary and because 
one was using phantom dust rather than dust. This is why the argument of the square root is classically forced to be 

non negative. Thus, classically wc can write H = C\Jl — ^ ah ^ Ch and know that the Taylor expansion of the square 
root 




(a) + 0([Q ab C a C b /C 2 } i ) (2.4) 



converges. One could now argue that square root and its Taylor expansion (|2.4[) provide classically equivalent starting 
points and therefore instead of defining the square root via the spectral theorem, one could quantise the Taylor 
expansion directly. To do the latter, one would need to order the terms in each order of the expansion symmetrically 
and one would need to define the operator corresponding to 1/C on a dense domain where the expansion is valid. To 
that end, notice the identity 

C e^Q+ c z + e z e^o+ 2 C + it C — it e— >o+ I 

where R(z) = (C + z)^ 1 denotes the resolvent of C at z $ spec(C). Thus if we can quantise (the smeared) density 
C(a) as a self - adjoint operator, then the resolvent R(±ie) is a bounded operator at finite e. The removal of the 
regulator will of course only be possible on a suitable dense domain of the Hilbert space, for instance on analytic 
vectors for C. 

In order not to deviate from the main thrust of the paper, we will neglect the higher order contributions in (|2.4I) 
by considering a subspace of the Hilbert space for which in the sense of expectation values < Q ab C a Cb ><C< C > 2 . 
We will come back to the higher order terms in a future publication. Thus for the remainder of the paper we will 
assume that 

#phy S JdaC(a) (2.6) 

s 

is a good approximation. For instance, for spacetimes close to a homogeneous one the approximation should be very 
good as one can see by expanding in terms of the number of spatial derivatives. Then the zeroth order contribution to 
the total diffeomorphism constraint vanishes identically so that -ff p hy S = V C 2 = C is true without any approximation 
in zeroth order. To define the sectors of validity of the approximation in the quantum theory is more delicate but one 
could use semiclassical perturbation theory [55]. See also [35] for related ideas. 
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2.3 Canonical transformation 



Using the gauge fixed formalism described above, where one has a preferred reference frame given by the Brown-Kuchar 
dust clocks, general relativity in terms of observables can be defined as follows: Let M. be an (infinite dimensional) 
symplectic manifold and (Q a b(o~), P ab (a)) a coordinate basis of A4. Q a b{o~) can be interpreted as a (physical, i.e. 
gauge independent) Ricmannian mctricLj on a 3-dimensional manifold E with local coordinate system {cr a }^ =1 defined 
through the readings of the dust fields and P ab (a) is related to the extrinsic curvature of E in a 4 - dimensional 
globally hyperbolic manifold M which in turn is isomorphic to R x E. E as a Ricmannian manifold carries geometric 
structures such as a unique torsionfree covariant derivative D or the Ricci scalar R(Q). Q a b{&) and P ab (a) form a 
canonical pair, i.e. the symplectic 2-form £ T*M x T*M is given by n := / dadP ab (a)AdQ ab {a). Further more let 

C 9 := - 7 f=[QabQod-\QacQbd}P ac P M --V^QR{Q) (2.7) 
^det Q I K 

G% := -2Q ab D c P bc (2.8) 

be C°° functions on A4, Xqb , X c g their respective Hamiltonian vector fields and k = 8irG/c 3 the gravitational coupling 
constant. These functions are the gravitational parts of the Hamiltonian and diffeomorphism constraint respectively. 
If we would work in the standard formalism of GR where q a b,P ab are the gauge dependent 3-metric and its canonical 
momentum then C 9 and C 9 (after the replacement Q ab — > q ab , P ab — ► p ab ) would generate gauge transformations and 
only quantities which are constants along the flow lines of their respective Hamiltonian vector fields would have a 
physical interpretation. However, in the enlarged phase space described in the last section, C 9 and C 9 are just the 
gravitational parts of the constraints and, without taking into account the parts coming from the dust fields, do not 
generate gauge transformations. Due to the gauge fixing and construction of manifestly gauge invariant quantities 
Qab(o~), P ab (a) we have a physical Hamiltonian that generates time evolution in the dust-frame, see equation ()2 . 2Q. 
For the purpose of what follows we introduce an orthonormal frame at each point a £ E and define the cotriad_| 
E l a through Q ab ( u ) — SijE^E^. It can easily be seen that this definition is invariant under $0(3) or more general 
SU(2) rotations, thus the cotriads can be interpreted as su(2)-valued one forms with k, ■ ■ ■ = 1, 2, 3 su(2)-indices. 
Furthermore define the triad E°- to be the inverse of the cotriad through EfE l b = 6 b , EfE{ = S{. 
Now we perform the following coordinate transformation: 

(Qab,P ab ) - (Ei,P?:=2P ab Ei) . (2.9) 

One can easily check that this transformation is indeed a canonical one, i.e. it leaves the symplectic 2-form invariant 
and thus (E^, P?) form a canonical pair 

{PZ(*),Et(a')} = K5 a b 5\5{a,a') . (2.10) 

However, the phase space spanned by (P^,E^) is larger than the one spanned by (P ab ,Q a b)'- Only the symmetrised 
product P ab = ^PjEj := j(P°-E b + P b Ej) survives when we invert the canonical transformation and want to 
transform back to ADM variables, which means that the canonical transformation (Q ab , P ab ) (E[, P" 1 ) is not one-to- 
one but one-to-many. To get rid of the additional degrees of freedom we must demand p] a Ef := |(P? E b - P b Ef) = 
which after some algebraic manipulations leads to a constraint 

<:.■ '^/'r^ . (2-H) 

the so called SU(2)-Gauss constraint!!. 

Modulo Gauss constraint the gravitational parts of the Hamiltonian and diffeomorphism constraints can be written 
as ° 9 = ZdifE) i E i E b S m ~ h E l E r\ P k P t ~ ^P-R{E) and C 9 = -E l a D b {e)P b respectively where R(E) denotes the 
Ricci scalar written in terms of the cotriad and D(E) is the covariant derivative compatible with the cotriad. 
So finally we arrive at a classically equivalent formulation of general relativity in the dust-frame, which we want to 
summarise in the following paragraph: 

7 This identification can be made precise by demanding appropriate boundary conditions for Q a b^P ab an <i specifying the function space 
in which they live. See 1181 for more details and appropriate boundary conditions. 

8 Again, we use a capital letter E l a to make clear that this is a gauge invariant version of the cotriad. It should not be mistaken for the 
densitised triad used in LQG. 

9 In analogy with the formulation of GR in terms of Ashtekar's variables we still call this the Gauss-constraint although it does not have 
the form known from Yang-Mills theories. 
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The phase space M! is spanned by an sit (2)-valued one-form E^ and its canonically conjugate momentum P?, a 
vector density which takes values in su(2) as well. The phase space is subject to the Gauss constraint 

Gi = e ijk P?E k a , (2.12) 

which means that only quantities which are invariant along the flow lines of its associated Hamiltonian vector field 
have a physical interpretation. In the gauge fixed formalism we have a physical Hamiltonian 

tf P hys = j da^J(Cs + O) 2 - Q ab (C a a + C™)(C§ + C fe m ) (2.13) 

that generates motion in dust-time and in our variables the gravitational parts of the Hamiltonian and diffeomorphism 
constraints are given by 

° 9 = KG lV? P j- -dct(E)R(E) (2.14) 

Hi 

CI = -ElD h {E)P\ , (2.15) 

As before C m and C™ denote all possible standard matter contributions that one might couple to gravity and we 
defined the de Witt-type supermetric G l ^ b := 4dc t( g ) [-^a-^b ^ — ^E l a E J b ]. Thus, in contrast to the formulation of 
general relativity in terms of Ashtekar's variables the Hamiltonian constraint looks similar to a standard Hamiltonian 
composed out of of a kinetic term oc P 2 and a potential det(E)R(E). 

The choice of these "hybrid variables" lies in between those of connection dynamics and geometrodynamics. Ge- 
ometrodynamics is not sufficient if we are interested in coupling fermionic matter. Furthermore, it is much harder to 
smear symmetric tensors over submanifolds in a spatially diffeomorphism covariant manner than it is possible for the 
one and two forms e, *P respectively. The drawback of using cotriads rather than connections, as already mentioned 
in the introduction, is that the holonomy transforms locally under SU(2) gauge transformations while the integral 
of the cotriad along one dimensional curves does not. However, in the algebraic setting that we choose in the next 
section, this will not pose any problem. 

For what follows we need the explicit expression of the curvature R(E)— term in the Hamiltonian constraint. The 

cotriad E l a together with the compatibility condition D a E b = uniquely defines a spin connection r* through 
D a El =: d a E l b + T c ab El + e^^r^E^ where T^ b are the components of the Levi-Civita connection. In terms of the 
cotriad and its inverse this spin connection can be written as 

Ti = -\e ijk E){El b -El a + ElE™E™ h ) , (2.16) 

where we used E l b as a short hand notation for d^E^. 

Using the above covariant derivative one can define a curvature tensor R^ ab in the obvious way aJ^l 

[D^D^Vj-ejuRtbVi (2.17) 

for every sit(2)-valued scalar function Vj on S. In terms of the spin connection this curvature tensor can be given 
explicitly as 

R l ab = 2d {a Y l b] +e lkj T k a T{ . (2.18) 

Finally the Ricci scalar in terms of cotriads can be written as 

R(E) - e jkl R l ab E«E b k = 2e jkl E?E b k d a T l b + (£^) 2 - E^E b k Y k a T{ . (2.19) 

In this form the Ricci scalar is a function of the cotriad E % a and its inverse Ef. Having a quantised version of this theory 
in mind, it will turn out to be more convenient to rewrite R in terms of the cotriad and its determinant only which 
can be achieved by employing the identity Ef = 2 det(-E) eabCe ijk^b^c- Thus the 3 terms in the preceding expression 



10 This definition is only meaningful if dim(E) = 3 where the fundamental and the adjoint representation of the gauge group have the 
same dimension. 



8 



take the form 

-EftfriT* = 4(det 1 (jg))2 c abc t def [KiEh - El J + EiKbl [Ef( E le ~ E k e J + E l a E l f>e ] (2.21) 
2e ijk EtE)d a V k h = {d J {E))2 ^^E\E) [E?^ + «d + 44,da] 

j abc def a'b'c' r?i jpi jpj rpk 

+ 2(det(E)) 3 e 3kl E c E f E b ,E d 

[E?{E l d>a E™ a , + E l eia E% td ) + E?E l ±a (EZ a , - E% <b )] . (2.22) 
This means that the "potential" -term in the Hamiltonian constraint is symbolically of the form 

/ / E 2 (dE) 2 E 3 (ddE) E 5 (dE) 2 

det ^) « w + w + (4)F ■ (2 - 23) 



In terms of the momenta P®, the cotriad E % a and its determinant the diffeomorphism constraint can be written as 
d = -ElD b P* = -EldtPP - e ijk EiT{Pl 



= -E l a d b Pt + \pt{El b - E{ a ) + —L-^e^e^E^P* [E l a (E l b>c E l cfi ) + E l b (E l atC - <J] .(2.24) 

Throughout this article we will work in units where c = 1 and the coordinates of E have the dimension of length i.e. 
[a] = cm. Thus, the metric Q a i, is dimcnsionlcss as well as the cotriads, [E l a ] = 1. The extrinsic curvature of E is 
measured in units of cm -1 which directly carries over to our canonical momenta, [Pf 1 ] = cm" 1 . The gravitational 
coupling constant carries units [k] = cm ■ kg" 1 and we will leave Planck's constant dimensionful, [H] = kg ■ cm. 

2.4 Coupling matter 

Now we want to couple matter to the gravitational sector. In principle one can allow all (standard model) matter but 
for illustrative purposes we want to restrict ourselves to the case of a minimally coupled scalar field of Klein-Gordon 
type in this article: The basic variables arc given by a scalar field $ and its canonically conjugate momentum II, both 
assumed to be in C°°(E): 

{n(cr),$(a)} = \6((T,a') (2.25) 

with the matter coupling constant A which we will set equal to one for the remainder of this article. 
In ADM-variables the Hamiltonian and diffeomorphism constraints in the matter sector are given by 



C* = nd a $ (2.27) 
respectively which can directly be rewritten in terms of our gravitational variables (E l a , P?) as 

c* = ;r^ n2 + ^^^^ (2-28) 

Ct = Ild a <i> . (2.29) 

In the above definitions we keep the factor of in the definition of the massive term to make sure that II, $ have the 
right dimensions: It can easily be seen that [$] = kg" 1 / 2 ■ cm" 1 ! 2 and [II] = kg 1 / 2 ■ cm" 3 ! 2 . 

The physical Hamiltonian, which generates evolution in dust-time for the whole system gravity plus matter is now 
given by 

iW = J dayj {C9 + O) 2 - Q ab (C 9 a + Ct)(C 9 b + Cf) . (2.30) 
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v + (0,0,1) 




v + (1,0,0) 

v- (0,0,1) 

Figure 1: A part of the graph 7 with vertices v,v ± (1,0,0), v ± (0, 1,0), v ± (0,0,1) and (outgoing) edges 
ci(v),c 2 (v),c 3 (v). 

After the canonical transformation in the gravitational sector the matter degrees of freedom only couple to the cotriads 
E % a and not to their momenta P 2 a . When turning to the quantum theory in section [3] we will choose a representation 
where the cotriads are quantised as multiplication operators as opposed to the LQG-representation where (smeared) 
densitised triads act as (non-commuting) right-invariant vector fields on SU(2). Albeit loosing the non-commutative 
structure (which might turn out to be important in the long run) this is one of the main advantages of working in this 
representation because having to deal only with gravitational multiplication operators in the analysis of the matter 
sector makes it much easier to analyse the QFT on CS limit in this framework. 



2.5 Discretisation 

Having rewritten general relativity in terms of the cotriad E 3 a and its canonical momentum P? we now want to explain 
how one can discretise this theory on a regular cubic lattice in 3 spatial dimensions. This will result in a classical 
dynamical system of an at most countable number of degrees of freedom and it is exactly this dynamical system which 
we will compare the quantum theory with that we will define in subsequent sections. 

Let 7 be a piecewise analytic graph, embedded in E and of cubic topology (that is, each vertex is 6-valent). For 
details concerning piecewise analytic graphs of the type usually employed in LQG see for instance [7j. For our 
purposes it is sufficient^ to think of 7 as a collection V(j) of points v G E, called vertices, together with analytic 
edges a[v] : (0, 1) — > E; s 1— > c^[v](s) such that cj[u](0) = v and c/[«](l) = v + 1. The requirement of cubic topology 
means that each vertex v has exactly 3 outgoing edges, that is, I,J,K= 1, 2, 3. Hence, each vertex v has exactly 6 
neighbouring vertices, which we will denote by v ± (1, 0, 0), v ± (0, 1, 0) and v ± (0, 0, 1), or simply as v ± / for the 
neighbouring vertex in direction /. 
Furthermore we will need the dual graph 7* whose faces are defined through a smooth embedding 

s a[v] ■ (0, 1) x (0, 1) -> E; si,s 2 i-> S CI [ v ](si,s 2 ) ■ (2.31) 

The embedding is chosen in such a way that each surface S Cl [ v ] £ 7* intersects 7 only in exactly one point p c , which 
is an interior point of S CJ [„] as well as c/ [v] . 

A discretised version of general relativity using the variables defined above should not depend on any background 
metric, since there simply is none. In order to define lattice variables that do not show any dependence on a metric 
one proceeds as follows: E 3 a is an su(2)-valued one-form, so we can naturally (that is without reference to any metric) 
integrate it along one-dimensional submanifolds of E, namely the edges. Pj 1 is an sit (2)-valued vector density, so we 
can naturally integrate it over the faces of the dual graph. Similiar reasoning holds for the matter variables: n is a 
scalar density, thus it can naturally be integrated over the volumina of the dual graph, and for $ as a scalar field of 

11 We will be rather informal here: Of course one should not consider paths as we do here but certain equivalence classes thereof. But in 
this article these details play a minor role so we will drop them for the sake of clarity. 
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Figure 2: A part of the graph 7 and the corresponding part of its dual graph 7*. 
density weight zero the natural thing to do is to just take the value of the field at each vertex. Thus, we define 



E){v) 


:= J dsdsc^i 









1 1 




f f 

:= / ds\ / dsa 







$0) 


:= <5>(a)\ a = v , 


U(v) 


:= J daU(a) 




V(v) 



dS c I [v] dS c I [v] 

dsi ds2 



£abcPi(S CI [ v ](si,S 2 )) , 



(2.32) 

(2.33) 

(2.34) 
(2.35) 



These discrete variables fulfil the Poisson algebrj^l 

{if («),£>')} 
{pi(v),P?(v')} = {E}(v),Ei(v')} 

{n(«),*(i/)} 
{n(«).n(i/)} = {#(«), *(t/)} 



, 
. 



(2.36) 



Poisson brackets between matter and geometry variables vanish obviously and the S vv i appearing in the above formulas 
is now the Kronecker delta. The Poisson brackets lose their distributional form exactly because we have chosen to smear 
the dynamical variables in dimensions which add up to dim(S) for canonical pairs. Due to the smearing procedure 
the discrete variables carry different units compared to the continuum variables, namely [-£}(«)] = cm, [-P/(t>)] = cm, 
[<f>(v)] = kg 1 ' 2 ■ cm- 1 ' 2 and [Il(v)} = kg 1 / 2 ■ cm 3 / 2 . 

Now we can define discrete analogues of the constraints from the last section by simply replacing 



Kir) 

$(a) 



E\{v) 

n(«) 



(2.37) 



and replacing the integrals J da by Riemann sums ■ Of course there is an ambiguity in defining these discrete ana- 

E V(a) 

logues just because there is no unique way of replacing the derivatives of E l a {cr): One could replace d a E l b {<j) — > VfE l j(v) 
or d a E l b {a) — ► \7J Ej(v) where Vf, VJ are the lattice forward and backward derivatives into direction I respectively 
defined through Vjf(v) := f(v + I) - f(v) and Wjf(v) := f(v) - f(v - I) for functions / : V( 7 ) -> T ( m our case 
either T = M. or T = sit (2) but one could have different topological spaces T). Both definitions (and others involving 



Here we fix the matter coupling constant A = 1. 
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not only next - neighbouring terms) are equivalent when taking the continuum limit. 
Remark: 

The algebra of the naively discretised Hamiltonian densities or any other quantities involving spatial derivatives that 
are replaced by difference operators will differ from their continuum limits even in the classical theory. This poses no 
consistency problem because these are no constraints anymore but true Hamiltonian densities, there are no anoma- 
lies. However, the densities will no longer be constants of the motion and one may want to improve the situation. 
If we would have constraints then one could fix the corresponding Lagrange multipliers by the method proposed in 
[471 EH S3 EH] • However, we are not in this situation and thus one would need to rather change the naive discretisation 
into a more subtle one as suggested in [20], possibly along the lines of [51]. Notice that we still do have a constraint, 
namely the Gauss constraint, but in our variables, in contrast to the connection variables, this does not involve spatial 
derivatives and thus there are no anomalies. 

3 Quantum theory 

In this section we want to define a quantum theory of General Relativity in the formulation developed in section [2] 
To be more specific, we will start with the algebra (|2.36| and explain how one can represent it on an appropriately 
chosen Hilbert space. First we will define a quantum theory on a purely algebraic level, that is, we will start with 
an elementary algebra of observables 21 labeled by an algebraic graph a and then construct a representation of this 
algebra on some Hilbert space Ti. defined on that graph. Thus, dynamics will be defined on an abstract graph, which 
does not know anything about embedding, differential or even topological structures of a manifold, since there is no 
manifold at this stage. In the framework of LQG many things do actually not depend on the embedding and this was 
one of the motivations in [14] to formulate an LQG inspired quantum theory purely on the algebraic level. However, 
if we want to interpret this theory as a quantum theory of General Relativity we need to make sure that there exists a 
large enough class of semiclassical states (and these states will carry information about the embedding of the graph a 
into the manifold S to be approximated) and it can indeed be shown that such semiclassical states exist (see section 
[4| . If the quantum system is in such a state the expectation values of geometrical observables will be sharply peaked 
on classical values and thus, we are getting classical discretised gravity as a classical limit of our algebraic theory. 

Now we want to define the quantum theory living on an algebraic graph and we will do this in a similar way as in 
[14j . First we need to define what we mean by an algebraic graph as opposed to the embedded, piecewise analytic 
graph 7 from section [2751 Let a be an oriented algebraic graph with N vertices, i.e. a countable set V(a) of N points 
("vertices") together with an N x N matrix a whose entries a vv ' , v,v' = 1, . . . , N take non-negative integer values. In 
general a VV ' is not a symmetric matrix for an oriented graph and the total number of "edges" connecting vertex v and 
vertex v' is given by c vv > := a vv > + a v ' V . To be more intuitive, one can think of a as all the information that remains if 
one starts with an embedded, piecewise analytic graph 7, consisting of N vertices and a vv > embedded edges between 
vertex v and v' and then drops everything that 7 knows about metric, differentiable or even topological structures of £. 
a neither contains any information about the braiding of its edges nor about smoothness or n-times differentiability 
usually associated to edges which are embedded in a manifold. 

For the remainder of this article we will be interested in the case N = H and graphs of cubic topology, i.e. each vertex 
v e V(a) is 6 valent and c vv > € {0, l}Wv,v'. Thus, each vertex v has three outgoing edges, which we will label by 
I,J,K = 1,2,3 for bookmarking purposes and three ingoing edges which will be counted as outgoing edges for the 
respective neighbouring vertex in direction —I, — J, —K. 

3.1 The algebra 

Now we define an abstract *— algebra 21 whose generators are labeled by vertices of the algebraic graph a and which 
implements the *— relations and canonical commutation relations that arise from the classical reality conditions (all 
discretised fields are real valued) and the Poisson brackets (|2.36|) : 

3.1.1 Gravitational sector 

Let a be the cubic algebraic graph with N = H as described above. With each vertex v of a associate a triple E}(v) 
of su(2)-valued quantities where I,J,K,... label the 3 outgoing edges at each vertex and i,j,k,... are sw(2)-indices. 
Furthermore associate a triple of su(2)-valued quantities (v) with each vertex where again capital letters denote 
outgoing edges and small letters are su(2)-indiccs. Let these basic variables be subject to the following algebraic 
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relations: 

[Ei(v),m(v')] = o 

[Pl(v),P?(v')] = 

[P/ (<;),£>')] = KSuSijS^ . (3.1) 
Furthermore we can define an involution on the algebra and will demand the (trivial) *— relations 

(E}(v)y =E\{v), {Pt(v)y=P*(v) . (3.2) 
Thus Ej(v) and P/(f) form an abstract *-algebra with the above relations which we want to denote as 2l 9 . 

3.1.2 Matter sector 

The same thing can be done for the matter degrees of freedom: With each vertex v of a associate a pair of real valued 
quantities $(w) and II(w). Let these variables be subject to the algebraic relations 

[$(v),$(v')] = 
[TL(v),U(v')] = 

[tt(v),$(v')] = 5 VV , , (3.3) 

and again define a (trivial) involution * through 

($(«))* = (n(w))* = IL(v) . (3.4) 

Again (f>(v) and H(v) together with the above relations form an abstract *-algcbra which we want to denote by SI 1 *. 
So the whole theory is described by an abstract * -algebra 21 := 2l 9 (g> 21^. 

3.2 The Hilbert space 

21 can be represented as an algebra of linear operators on an ITP Hilbert space := ® V TL V where H v = Ti.% ® H™. 
where the individual Hilbert spaces on each vertex, "H% and H$, are chosen as explained below. Note that this ITP is 
different from 

H%®Hi:={® v Hl)®{® v Ht) , (3.5) 

because of the non associativity of the infinite tensor product. The bracketing chosen by us is preferred because it 
allows to construct gauge invariant quantities. More information about the ITP construction can be found in appendix 

m 

3.2.1 Gravitational sector 

Let 1-L% := Z/2(K 9 , d\i) be a Hilbert space associated to each vertex v with Lebesgue measure \i on R 9 and consider the 
ITP Hilbert space H® ■ On the algebra 21 s can be represented by choosing a representation p : 21 s — > C(H®) in 
which E}(v) acts as a multiplication operator and P^ (v) as a derivative operator, 

E](v) ■ HE) := [p(E}(v))i>] (E) := E}(v) ■ ip(E) (3.6) 

P!(v)-iiE):=[p(P!(v))i>](E) := U%^^^(E) , (3.7) 

for all il>(E) € 7i|, and lp — \/kK is the Planck lengtlS Note that the derivative in the second line is indeed just the 
partial derivative and not a functional derivative. As operators on 7i|, the cotriad E}(v) and its conjugate momentum 
P/(u) fulfil the commutator relations 

[P^v),Pj(v')] = iPpSuSySn, . (3.8) 

In what follows we will mostly omit the hat above the operators and use the same symbols E\(v), Pf {v) for elements 
of the abstract * -algebra 2l 9 and their representatives as linear operators on the Hilbert space Ti% but it should be 
clear from the context which one we mean. 



J We need to include a factor of l 2 p into the definition of the operator P?(v) in order to get the same dimensions as its lattice-analogue. 
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3.2.2 Matter sector 

Again, a natural representation of 21* is given by the ITP representation on 7i* := ® v 7i% where each vertex labels a 
Hilbert space Tif := dv) with v being the Lebesgue measure on R. We choose a representation such that 

:= <Z>(v)tp(v) , (3.9) 
E(«)-V(*) := ihd Hv} ^) , (3.10) 

for ip £ TCL . As linear operators on Tt^ these fulfil the commutator relations 

[%),$W]=iM„„- . (3.11) 

3.3 The Hamiltonian Density 

In section \2.2\ we saw that, given that — ^ Ch is small, one can approximate the physical Hamiltonian by flphys ~ 

/ da (C 9 (a) + C*(cr)) plus higher order terms. So in order to define physical dynamics in the quantum sector, the first 
s 

step is to have a well defined quantum operator C associated to the non dust contributions of the classical Hamiltonian 
constraint J da (C 9 {a) + C*(er)). For higher order corrections we will need the diffeomorphism constraint as well, so 

s 

for completeness we will also give a quantum operator C a (v) associated to the non dust contributions to the classical 
diffeomorphism constraint density. 

The crucial step in defining these operators is to be able to define operator analogues for functions involving inverse 
powers of det(-E). In this section we will just formally define these operators, in the next section we will show in detail 
that they are in fact symmetric operators on a dense subspace of W<g>. 

All the operators will be of the form O = O 9 <E) O*, where the first part, O 9 , acts only on Ti% and the second part, 
O*, acts only on H®. 

3.3.1 Gravitational sector 

As the operator analogue for the (integrated) gravitational part of the Hamiltonian constraint we defincF^I 

& ■■= (3-12) 
Cpot := Cpoti + Cpot 2 + ^pot 3 (3.13) 

<%n ■= «EE[^!W]^K 9 4W] ' (3 ' 14) 

V I.J 

^ := A K detE(v) £ e IJK e LMN E k (v)E l L (v)(W+E k (v))(W%E l M (v)) , (3.15) 

^ ' I,J,K,L,M,N 



IJK LMN 



r 9 ■= V 

pot2 ' 4ndctE(v) ^ 

y ' I,J,K,L,M,N 

[E k K (v)(VjEUv)-VtE k j(v))+EUv){V+E k K (v))] x 

x [4(«)(ViEf(») - VjE k M (v)) +E\{v)(Vl I E l N {v))} , (3.16) 

v ' I,J,K,L,M,N 

x [E l M (v)VjViE l j(v) + (V+E l j(v))(V+E l M (v)) + E l ,(v){VJV + L E l M {v))] 



UK LAIN OPQ 



■ j uE i K (v)E i N (v)E 3 P (v)E%(v) x 



2n(dctE(v)) i <JiK ,l'm,n,o,p,Q 
x [EJ(v)[(V+E l L (v))(V+EZ(v)) + (VjE l M (v)){V+E%(v))] 
+ EZ{v){VjE l L {v)) [(y+EJ(v)) (V+Eg»)]] . (3.17) 



14 For the most part of what follows we will again drop the hats over the constraint operators when there is no ambiguity in mistaking 
them for the classical functions. 



14 



id^ is the identity operator on . The multiplication operator in the kinetic term is given by 



(3.18) 



and by and V 7 we denote the lattice forward and backward derivatives into direction I respectively defined 
through Vff(v) := f(v + I) - f(v) and VJf(v) := f(v) - f(v - I) for functions / : V(a) -> su(2). 

The operator Cf(v) = Cf(v) ® id,/,, associated to the classical gravitational diffeomorphism constraint density, can be 
defined as 



1 



(3. 



J,3 



■I.J 



- E 

8 ^ 

JKLM 

ijkl 



KLM 

e eijk 



L4(,)4( B )[4(,,)((V+^(«)) - (V+E l K (v))) + E l J (v)((V+E l I (v)) (VjE l K (v)) 



det(E) 



det(E) 



Pf(v) 



Here we used the notation V~[Pf(v) := Pj (v + I) — Pj 7 (v) = iipd^ („ + j) — ^p^E j (v)- This operator, besides the 

terms proportional to V^Ej(v), causes Cf(v) to act not only on vertex v but also on its neighbouring ones. 
Finally the gravitational 5f/(2)-Gauss constraint operator can be defined as 



G[A] = id ®i Yl €ijk[Pttv)EHv)Wv) + EHv)Mv)li(v) 



(3.20) 



veV(a) 



for A : V(a) — > SU(2); v h-> A,-(u). 



3.3.2 Matter sector 

As operator analogues for the (integrated) matter Hamiltonian constraint and the matter diffeomorphism constraint 
density we define: 



E 



I, J, K, L.M.N 



2 det 



e /^ e ™^ (w)£; j. {v)E k M{v)E i N{v) g, (V+$(«))(V+*(«)) 



+ det£(u) ® ^5-$ 2 (v) 
W ff ® y [(Vt*(»))«*(«) +9* w (V+$(«))] 



(3.21) 
(3.22) 



One can rewrite (|3.21[) in a slightly different way which will become more convenient for the analysis of the matter 
sector which we want to perform later on 



c* = - Y Wv) - 1 . M (v) 

o ^ y ' det£(u) y ' 



vGV(a) 

with the lattice Laplace-Bcltrami operator A defined as 



(3.23) 



A 



1 



det E(v) 1 j 



E v 7 



K,L,M,N 



1 



JKLJMN j^k 



2 det E{v) 



E" K (v)E L (v)E M (v)E' N (v)V: 



(3.24) 
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Here we used that V|" = (V7) if interpreted as an operator on l 2 (V(a)), the Hilbert space of square-summable 
functions on the set of vertices V{a). This Hilbert space, l2(V(a)), should not be mistaken for the ITP Hilbert space 
which we based our quantum theory on. But it will play a role as the "one particle" Hilbert space in the construction 
of an appropriate Fock space. 



3.4 Domains of Definition 



In order to show that the operators written down in the last section are actually well defined symmetric operators on 
a dense subspace of H.® we proceed as follows: As can easily be seen the parts acting on will cause no problems. 
For each vertex v we can define the matter part of (|3.21[) and (|3.22p on the Schwarz space S(M) of smooth functions 
of rapid decrease which is a dense subspace of Tt% = L 2 (R 5 dv). 

In order to show that the operators (|3.12j) , (|3.19|) as well as the gravitational part of (|3.21|) and (|3.22p is well defined 
it is sufficient to show that A n (v) :— (ggtpSTOJT" ^ or n = 1j 2 can be defined on a dense subspace of H.% = L 2 (R 9 , dfj,) 
and that the entire operators are symmetric on that subspace. 

Let us analyse the determinant of the cotriad in more detail: The operator E\(v) (associated to the cotriad) is densely 
defined on S(R 9 ). so we need to show that this holds also for powers of the inverse of its determinant . If we view 
the determinant det : R 9 — > R; E\(v) t—> dct(E)(v) as a real-valued function on R 9 then the singularity structure of 
dct(-E)(t>) can mos t easily be analysed if we perform a singular value decomposition as follows: 

E =: LDR , (3.25) 

where L, D e SO(3) and D is a diagonal 3 x 3-matrix with real eigenvalues^!. Explicitly these matrices are given by 



'E\ E\ 



E 



L 



El 
El 



cos (/3) cos (7) 
sin (a) sin (f3) cos (7) + cos (a) sin (7) 
- cos (a) sin ((3) cos (7) + sin (a) sin (7) 



— cos ((3) sin (7) 
- sin (a) sin ((3) sin (7) + cos (a) cos (7) 
cos (a) sin ((3) sin (7) + sin (a) cos (7) 



sin ((3) 
— sin (a) cos ((3) 
cos (a) cos ((3) / 



(3.26) 



D := 




— cos (9) sin (ip) 
- sin (S) sin (9) sin (ip) + cos (<5) cos (<p) 
cos (5) sin (9) sin (ip) + sin (5) cos (ip) 



\ 



sin (6) 
— sin (5) cos (9) 
cos (5) cos (9) J 



( cos (9) cos ((p) 

R := sin (5) sin (6) cos (<p) + cos (5) sin (<p) 
y — cos (S) sin (9) cos (ip) + sin (6) sin (<p) 
That means we perform a coordinate transformation 

p : M 9 -> R 9 ; (E$, Ef, Ef, E%, E%, E|, E\, E%,E$) h-> (A x , A 2 , A 3 , a, /3, 7, 5, 9, 4>), where the range of the new coordi- 
nates is a, "f,S,(f> £ (0, 27r], [3,9 £ (0, 7r] and Ai, A2, A3 G (—00, +00) and the Jacobian of this coordinate transforma- 
tion can be computed to be 



det($) := det(<9</Vd(A, a, .., 4>)) = cos(/3) cos(9)(~\ 4 3 \l + \fX 2 2 - X^Xj + X 2 3 X^ + X 4 3 Xf - XjX*) 
L, R G SO (3), thus the determinant in this coordinate system is simply given by 

det(E) = X 1 X 2 X 3 . 
So if we compute ||^4 rl (w)i/'||L 2 for ip £ 5(R 9 ) we realise 

Ai) 



\A n (vM] 



dfl |cos(/3)cos(6>)| dXidX 2 dX 3 



Af(A|-Al) + A|(Al-Af) + Af(Af 



\2n \2n \2n 
A t A 2 A 3 



*(«,A 2 



(3.27) 



(3.28) 



(3.29) 



where we used £1 as an abbreviation for all angles and dfi := dadf3d / ydSd9d(/). The denominator in the last integral 
is zero if either Ai, A2 or A3 is zero whereas the numerator has its roots at Ai = ±A 2 , A2 = ±A3 or Ai = ±A3. This 



15 We could have absorbed the signs of the eigenvalues into L, R by having them take values in 0(3). But it is more convenient to have 
them take values on the whole real axis. Notice that we do not insist on definite sign of det(E) on the whole phase space. 
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means that the argument of the integral (including the factor coming from the measure) diverges for A^ — ► and thus 
A n (v) := ^ Act ^)(v)) n d° es no * nave finite norm on all ip G S(M. 9 ). 

To cure this problem we need to change the domain of the operator A n (v), and the heuristic idea is to choose as 
its domain functions ip G <S(R 9 ) with the additional condition that ip vanishes faster than any power of A^ in a 
neighbourhood of the coordinate hypersurfaces defined through Aj = 0. 

In order to do this we will modify the function ip in a way that they have the required properties and still approximate 
the original functions well. For this purpose we adopt a standard regularisation procedure. Let us introduce a 
characteristic function xv with V being the interval [— 2e, 2e\. \v is equal to one if x lies in V and zero otherwise and 
is not smooth. Now consider the family of functions 



Pe(x) 




(3.30) 



The constant C e is determined by the condition J m dxp e = 1 and is explicitly given by C e 1 = e J exp(— 1/(1 — y ))dy. 

This function is smooth, has compact support since supp p t C [— e, e] and on all of R we have p e > 0. We use p £ to 
regularise the characteristic function xv and obtain a family of smooth versions of the characteristic function which 
is equal to one if x lies within the interval [— e, e] and vanishes for all |x| > 3e. Let us denote the regularised function 
as Xe then 



2 c 



Xe(x) 



dyxv{y)pe(x - y) 



dype(x-y) 



(3.31) 



With this definition of Xe we can easily show that it vanishes if |x| > 3e. Since the aim is to set the functions at the 
singularities to zero in a smooth manner we define the function 







for \x\ < e 



+2e 



s e (x) := 1-Xe(x)= { 1- J dyp e (x-y) for e < \x\ < 3e 

-2e 

1 for Id > 3e 



(3.32) 



Notice that < Xe 5 s € < 1. This function is smooth by construction and s e (0) = 0. Thus, we can use this function to 
modify our original functions ip G 5(M 9 ) such that they vanish at A.; = in a smooth way. Define 



0,Aj h-> ^£(^,A 4 ) = s e (Ai)s e (A 2 )s e (A 2 )V'(fi ) Ai) 



(3.33) 



for any ip £ 



The product of two smooth functions is still smooth, hence the set T := {ip e } is a subset of S(l 



and one can easily see that T is dense in <S(R 9 ) and therefore also in 



in the L 2 -norm: 



||^-^|h 



< 



+00 



+oo 



dn J dXi I d\ 2 I d\ 3 dct^(n,x t ) - ^ e (n,\i)\' 

— oo 
+3c 



-oo 
h3e 



h3e 



dfl J dX 1 J dX 2 y dA 3 det A,)| 2 |l - Se (A 1 )s e (A 2 ) Se (A 3 )| : 

-3e -3e -3c 

1 

+3c +3e +3e 

dfl [ dAi I dX 2 J dA 3 dct$|i/'(^,A l )| : 



-3e -3c 



< l max V¥(? / dn\cOs((3)cOs(d)\) (X 



(3.34) 



where I max is the supremum of 



dct * 



|^| in the compact set |A»| < 3e which is assured to be finite for every 



I cos(/3) cos(0) I 

smooth function of rapid decrease. Thus, we can approximate every ip G S through some ip e G T to arbitrary precision. 
Hence, A n with domain T is a densely defined operator in i 2 (R 9 ) and therefore all the operators defined in the last 
section are actually well defined and even more they leave the domain invariant. 
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4 Coherent states for the gravitational sector 



In this section we want to construct coherent states for the gravitational sector: The gravitational Hilbert space 
introduced in section EOl is at each vertex of the algebraic graph given by TL 9 V = L 2 (M. 9 , dfj.) and thus consists of nine 
copies of usual Schrodinger representation Hilbert spaces. We will construct coherent states for each vertex Hilbert 
space H.Q and the coherent state for the complete gravitational sector are then given by the infinite tensor product of 
the vertex coherent states. Since the quantum theory is formulated on a purely algebraic level, the coherent states will 
also be the bridge between the algebraic quantum theory and the embedded classical theory of General Relativity. 

In order to construct coherent states we first need to choose an embedding X of our algebraic graph a in to a 
given manifold £ and we call its image 7 := X(a). We will choose embeddings X such that 7 is dual to a certain 
triangulation 7* . Thus, for each embedded edge X(c) there is a face S c in 7* which intersects X(c) only in an interior 
point p c of both S c and X(c). Next we need to choose a classical cotriad E (suppressing indices) and an su{2) - 
valued vector density P . With the manifold, the embedding and the classical data (E , -fo) we can define the smeared 
quantities 

E\ Q {v) = J E 

Pl (v) = J e abc da a Ada'P^a) (4.1) 

We follow the complexifier method for coherent states introduced in [12"1 02 [lOj [IT] . We complexify the configuration 
space R 9 to C 9 by introducing 

Z>):= -1= iPf Q {vj) (4.2) 



From now on we will suppress indices and denote the classical data simply by Z(v) = (E (v), P (v)). This complexi- 
fication can be obtained from a classical complexifier of the form 

with 

00 

z W = J2^{ E3 M c hn) (4-4) 



where {., .}(„) denotes the iterative Poisson bracket defined as {Ej(v), /}(o) = E J j(v) and then iteratively by 
{Ej(v), /}( n +i) = {-Ej(u), {£j(u), /}(„)}■ Similarly we can define a corresponding complexifer for the quantum 
configuration space 

C := Y (i£^—)(ul—^—]d jk S JK (4.5) 

from which we obtain the form of the annihilation and creation operators 

that satisfy the algebra [Aj(v), (A^Y (v 1 )] = 8 vv i8' jk 8jK- Note that the characteristic length that enters into the 
definition of the coherent states is £ p in our case and not as for the harmonic oscillator I = ^Jfijmu. The reason 
for this is simply that the squared Planck length £ p also occurs in the definition of the momentum operator. One of 
the defining properties of coherent states is that they are eigenstatcs of the annihilation operator AK with eigenvalue 
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Z , the classical phase space point that they are labeled by and around which they are sharply peaked. They can be 
expressed as 



|l&z) := ^exp(--|Z| 2 )— [(ASM) ] |0)=exp(--|Z| 2 jexp(z(^)n^))|0) (4.8) 



where |0) is the vacuum state defined by Aj(v)\0) = 0. Further properties of coherent states are that they satisfy a 
resolution of identity, they form an overcomplete basis, and that they are not mutually orthogonal but (ipz 1 , ipz) = 
exp(-|Z-Zf). 

In position space representation we get using Z(v) = (Eq(v), Pq{v)) 



M&M) = ^== cxp (-^Po(v)E (v) + ^P„(«)2#(v)) cxp (-± (e^v) - E Q {v)f^ (4. 



9) 



An easy calculation shows that 4>z(Ej) is indeed an eigenstate of A 3 j(v) with eigenvalue Z(v). The states ipz(Ej) 
are the ordinary harmonic oscillator coherent states and thus we can construct the coherent states associated to each 
vertex as a product of nine harmonic oscillator coherent states and the total states for the gravitational sector then 
as infinite product of states over all vertices of the graph 

3 

Mv) ■= II ME j j(v)) € U% and := (® v ^z(v)) G H% (4.10) 
j,.j=i 

The fact that -0Z is an eigenstate of the annihilation operator A J j with eigenvalue Z has the consequence that it 
approximates expectation values of Ej and Pj 7 semiclassically well. Precisely, we have 

(i>z , E j j $ z ) = E , (i> z , Pf i>z) = Po (4.11) 

However, when we consider the inverse operator 1/ det(£'j) it turns out its expectation value with respect to ^ z 
diverges. A simple way to see why this happens is to take a cotriad configuration at one vertex where all but two 
triads are vanishing. Then we have at least one integral of the symbolic form J R e~ x /(x — Xo)dx which does not 
converge. For this reason we have to introduce modified coherent states along the lines of the modified Schwartz 
functions of rapid decrease in section f3. 41 Thus we define the modified coherent states by means of the function s e in 
equation (|3.32[) as 

3 

r Z (v) := II s e(det(£>))V^(i?>)) and ^1 := (®v^|(v)) (4.12) 

More precisely 

?M(w):M 9 ^M; SI, A, ^ ^(fi, \A = Se (A 1 )s e (A 2 ) Se (A 2 )V'z(r!, A,) (4.13) 
Using the same steps as in equation 13.341 we obtain 

Wz-i>z\\ < l max V¥?(^J dSl\cos(0)cos(e)\y ocei (4.14) 



is the supremum of 



dct* 



\ipz\ in |Ai| < 3e which is finite for smooth functions. Hence, 



s(/3) cos(6») 

for sufficiently small e we can approximate the states ipz by the states ij) z for all values of the classical phase space 
point label Z . Furthermore, because the coherent states have the property of being sharply peaked around Z : the 
absolute values of ij)z that appears in I max in the estimate above becomes smaller and smaller the farther the classical 
phase space point Z is away from a singular configuration. Consequently as long as we are mainly interested in states 
which are not peaked around singularities the approximation works very well. Indeed we are actually interested in 
those states because we want to use these states to mimic a classical geometry background when investigating the 
semiclassical limit of this theory understood as QFT on curved spacetimes in section [5j 
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We will now verify that the expectation values for Ej and Pj 7 are also well approximated for sufficiently small e. 
The modified states are no longer eigenstates of A 3 j with eigenvalue Z but we get a correction involving the derivative 
of s e and the original state ipz- 



Wz + 



£ p ds e (det(E 3 j) 
dE\ 



1p2 



(4.15) 



This correction term only appears when | det(Ej)\ < e since otherwise s e (det(Ej)) = 1. For the error of the expectation 
value of Ej we therefore get 



(a e (det(Ej)) + l) $z(E j j)E j j (a e (det(Ej)) - l) M&j] 



dfid 3 Adet($) ((s e (A 1 ) S£ (A 2 ) Se (A 2 )) 2 - l) ^(A*, fi)£^(Ai, fi^A,-, fi) 



/ dfi J d\i J d\ 2 I dA 3 |det($)||(( S£ (A 1 ) Se (A 2 ) Se (A 2 )) 2 -l)||v;z(A l ,r!)^(A l ,r!)Vz(A l ,r!)| 

[S 3 /Z 2 ]x[S 3 /Z 2 ] -3c -3e -3e 
3e 3e 3e 

< J dQ J dXi J d\ 2 J dA 3 |det($)||^(A i ,n)^(A i ,f2)^(A i ,Q)| 

[S 3 /Z 2 ]x[S 3 /Z 2 ] -3c -3e -3e 



,3,3 



dfi |cos(/3)cos(0)| oc e 



(4.16) 



with X max is the supremum of 



dct* 



s(f3) cos(9) 



\ l 4>zE J A\i)ifjz\ in |Aj| < 3e. We used that E, s e commute and in step three 



we performed again a transformation to the (A, fi) coordinates and considered Ej as a function of these variables. 
Similarly we obtain for the momentum operator 



(r z ,p-rz)-^z,pf^z) 



< 



([VI - M A J V4> + Wz - i&z]) 

([V>l ~ M + s e [PfiPz]}) + (P?tPz, Wz ~ M) 

(<PZ, \[[P]{Se - I) 2 ] i>Z + [S e - 1] [Pjtyz]]) + (P/^Z, [* 

dfid 3 A|det($)| (s e (AiK(A 2 )s e (A 2 ) - l) 2 lif 



(4.17) 



1] i\>z t 
IVz(Ai,fi)| 2 



dnd 3 X |det($)| |Vz(Ai,f2)| 



Ai.fi) 



^z(Ai,n)| s e (A 1 )s e (A 2 )s e (A 2 )(4.18) 



where we introduced the abbreviation dEj(Xi, fi) = (dXi/dEj)dXi + {dSli/dEj)dili for the transformed partial deriva- 
tive. Here in the last step we have performed an integration by parts. As before, both terms above are only non 
- vanishing in the compact interval A, G [— 3e,3e]. Restricting the integration domain and using the upper bounds 
|s e |, |1 — s e \ < 1 we see that both terms are bounded by [6e] 3 [J dCl | cos(/3) cos(0)|] times the supremum of the respec- 
tive integrand on the domain |Aj| < 3e. Since that supremum does not increase if we lower e we see that the bound 
can be arbitrarily small. Thus also for the expectation value of Pj 7 we are arbitrarily close to the original value when 
using the modified coherent states. Finally, let us compute the expectation value for the operator A n := [1/ det(Ej)] n 
with respect to the modified states. To display the exact formulas requires a lot of notation, so we satisfy ourselves 
with giving the flavour of it by considering a simpler example which nevertheless contains all the essential features. 
The exact calculation proceeds completely analogously. We model the expectation value of A n by 



< 



>:= 



dx 

7^ 



-(x-x ) 2 /5< 



(4.19) 
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and want to compare its value with l/|xo| n . We write 



< -p; >:= h + h := 
and have due to \s e \ < 1 



e<|a:|<3€ 



dx 



\x\>3e 



dx 

7^ 



-(x-x ) 2 /8 2 1 



I'll < 



< 



< 



1 



dx 
3£ dx 



-(x-x a f/8 2 _ 



e -(x-\x \f/S 2 +e -(x+\x \) 2 5 2 



< e -(\x a \~3e) 2 /S 2 +e -(\ X0 \+e) 2 /S 2 



where we have assumed that |xo| 2> e. Next 
1 



h- 



:= -h +h~- 



dx e _ (x _ xo) 2 /5 2 _1 



\x\<Zt 



\x \ n 



|a;|>3e 



dx 



\x\ n |ao| T 



where in the last integral we have exploited the invariance under x — > —x. We can estimate 



3e 



To estimate h consider 



s 2 



n ln(|x|) 



(4.21) 



(4.22) 



(4.23) 



(4.24) 



The extremum y of / satisfies 2x(x — \xq\) = —nS 2 i.e. y± = [\xq\ ± \/x\ — n5 2 }/2 assuming |xo| 3> 6. Hence 
y + m \xo\, y- w n5 2 /4|x |. Since /"(x) = ~2/5 2 + n/x 2 wc find 

f»(y+) « -2/£ 2 + n/x 2 < 0, /"(*,_) « l/-5 4 [-25 2 + lG^/n] > (4.25) 

Thus only y + is a maximum and it also satisfies \y+\ > 3e. Accordingly, a saddle point estimate of I4 yields 

r °° dx 



00 dx 



< 



\x \ n 



^5 
\xo[ 

y+ 



J(y+)+f"(y+)[x-y+] 2 



7M T 



-(x-\xo\V/S- 



/\xo\ n 



\xo\ 



- 1 



- 1 



1 n 2 S 2 



\x \ n 4|i | 



Summarising 



< 



1 



> \x \ n - 1 



<(\h\ + \h\ + \h\)\x \ n < 



n 2 S 2 



(k 



-\x \ 2 /5 2 



Mx Q \ 2 V?<* 



4|sq|" 



-\x \ 2 /S 2 



(4.26) 



(4.27) 



In our application 6 = £p is fixed while e — > 0. We have already seen that |xo| ^ e,<5 in order that (|4.27|) is a good 
estimate. Let y = \xq\/5 then the third term in (|4.27|) is of the form 



ra-l 



y n e-y 2 ' 2 



(4.28) 



which should be small as compared to the fluctuation term n 2 5 2 / A\xq\ 2 . This shows that we cannot let e — > for our 
approximate coherent states. However, choosing e.g. e — d we see that for n of the order unity already for y of the 
order y/n the third term in (|4.27[) is subdominant. 
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5 Analysis of the system 



5.1 Born— Oppenheimer approximation 

To analyse the interplay between gravitational and matter dynamics and to understand how a (matter-) quantum field 
theory on a fixed classical spacetime can emerge from the fully quantum gravitational setup we want to employ an 
approximation scheme of Born-Oppenheimer type. As we will see in the following for such an approximation scheme to 
work it is essential that the geometrical variables are quantised as multiplication operators, not as derivative operators. 
This is the main motivation for not using the LQG-representation but the one that we constructed in previous sections. 
This representation captures some of the features of LQG (or the recently proposed generalisation coined Algebraic 
Quantum Gravity [331 Q31 US] ) , especially the fundamentally discrete graph-like structures which appear in the very 
definition of the kinematical Hilbert space. In order to be able to use a Born-Oppenhcimcr approximation scheme we 
have chosen a representation of the gravitational Poisson algebra which is commutative in the geometrical variables, 
unlike the non-commutative flux operators in ordinary LQG. Thus, the theory constructed in this paper is a close 
relative to LQG when defined on an algebraic graph which shares some of its features (based on graph-like structures) 
but deviates from it in other essential characteristics (abelian geometrical operators). Of course, the long term goal will 
be to generalise the methods derived in this section to full LQG also taking into account the non-commutative structure 
of the holonomy-flux algebra. But, as we will see below, due to the specific assumptions made in the derivation of 
the Born-Oppenheimer approximation such a generalisation is not straightforward and will require substantially more 
work. One step into that direction is the construction of a non-commutative flux representation for the holonomy-flux 
algebra, where the flux-operators act as multiplication operators on an appropriately defined Hilbert space and the 
non-abelian nature of the flux-operators is taken into account through a specific ^-product on that space |52| . 
The use of Born-Oppenhcimcr approximation schemes in the context of quantum gravity has a long tradition, especially 
in the older Whcclcr-deWitt approach (for a historical account see for example Kicfer's book [27] and references therein, 
for a recent application in the framework of spinfoam models see [28]). Most of the results so far were obtained on a 
rather formal level, because the Hilbert space on which Wheeler-deWitt theory is assumed to be defined is not known 
at all, but the application to minisuperspace models for FRW-universes gives reasonable results and can shed some 
light on conceptual questions such as the often discussed problem of time. 

Concerning this issue we are in a better position: The Hilbert space of the quantum gravity theory we consider is 
explicitly known and all the operators we use are densely defined symmetric operators on that space. Therefore the 
derivations which follow are all well defined and work beyond the formal level - at least in principle when neglecting 
practical difficulties in actually performing the calculations. Moreover, we are advocating a reduced phase space 
quantisation of General Relativity where the preferred reference frame is given by Brown-Kuchar dust fields. Thus we 
have a true physical Hamiltonian that generates evolution in dust-time, not only a vanishing Hamiltonian constraint. 
The problem of time is therefore solved already on the classical level and does not import any extra obstructions into 
the quantum theory. 



5.1.1 General Framework 

For simplicity we will first explain the Born-Oppenheimer approximation using a very simple quantum mechanical 
system and later generalise to the quantum gravitational setup. Assume that our classical phase space is four di- 
mensional and coordinatised by two pairs of canonically conjugate variables (Q,P) and (q,p) (a generalisation to an 
arbitrary number of configuration variables is straightforward but will not be considered here for pedagogical reasons). 
Further, assume that dynamics is generated by a Hamiltonian of the form 

H--=Sf + ^ + V(Q, q ) , (5.1) 

where M and m are some parameters describing the system under investigation (in the case of two particles these will 
be their respective masses) and V(Q, q) is a potential term which depends only on the configuration variables, not on 
their momenta. 

This system can be quantised on a Hilbert space H := TiP Cg> Tt q where HP and TC q are taken to be the usual spaces 
of square integrable functions over the real line. On Tt one defines operators in the usual way as (q^)(Q, q) := 
q^{Q, q), (Qtf )(Q, q) := Qtf (Q, q), (ptf)(Q, g) := ihd g *(Q, q), (P*)(Q, q) := ifrd Q ^>{Q, q) for functions *eH. Thus, 
the Hamilton operator of the total system is given by 

ds =-^-^+m?) , (5.2) 



22 



where Aq, A q denotes the Laplacian with respect to Q and q respectively. 

In general, depending on the exact expression for the potential term V(Q,q), it can be difficult to obtain solutions to 
the full eigenvalue problem 

H^ a =A a q a , (5.3) 

where a is a set of quantum numbers labelling the eigenvalues A Q and the corresponding cigcnfunctions However, 
under certain assumptions it is possible to obtain approximate solutions: Assume that (i) m <C M, which means that 
there are two clearly separated energy scales in the problem at hand, and (ii) we have some control over one half of 
the problem, namely we can solve the eigenvalue problem for the light variables 

H{Q)Xi{d\Q) = HQ)Xi{d\Q) , (5.4) 

where H{Q) := — -l^A q + V(Q,q) is an operator acting on the space Ti q that depends on an external parameter Q 
and Xiil't Q) (labeled by a set of quantum numbers i) are assumed to be square integrable in H. q and also depending 
on Q in a parametrical way. So in a sense this means we can solve 15.41 for each Q separately. 

Let us first stick to the first half of the problem and analyse a little further equation l]5.4[) : Assume that we can 
solve the eigenvalue problem for each value of Q and we have obtained eigenfunctions x(g; Q) £ TL q and eigenvalues 
A(Q) that depend on the external parameter Q. Then a priori, it is not clear how \ arid A behave if we change U L 
Berry [S3] showed that this question can be answered in a nice geometric way: Under the assumption of adiabaticitM 17 ! 
(i.e. assuming that the change in Q is so slow that the quantum system at time t is in a state x(<liQ(t))) it turns 
out that when moving the quantum system around a closed curve C in parameter space (in our simple example the 
parameter space would just be R but the same phenomenon holds also for higher dimensional parameter spaces) the 
wave function x picks up a non-trivial phase factor. That means when the system is prepared to be in an energy 
eigenstate Xi(0) for t = 0, after having followed the loop C in a time interval T it turns out to be in the state 

T 

Xi(0). The second exponential is just the usual dynamical phase factor 



i>(T) -exp [i 7j (C)]exp 



IfdtXiiQit)) 



o 



but the first one, explicitly given by 

li {C) = inj>{ Xi {Q)\d Q \xi{Q)) q , (5.5) 
c 

is indeed non-trivial. This phase factor, known as Berry's phase in the literature was experimentally confirmed in 
many experiments and can be used to explain a variety of physical phenomena, including for example the Aharonov- 
Bohm effect (see [54] for more information on this broad topic). From a geometrical point of view the situation 
is quite interesting (see for example |55j and references in there for a modern explanation using the language of 
fiber bundles), because it turns out that A :— iH(xi \ 9q \ Xi) q transforms as a t/(l)-connection under coordinate 
transformations in parameter space. Generalised to a higher-dimensional parameter space this connection is given by 
:= ih(xi I I Xij m a local coordinate system Q M . Thus, exp [^7i(C)] is nothing else than the holonomy of 
a [/(l)-connection along the path C in parameter space. 

One essential point in the derivation of the above formulas is that Q is treated as an external classical parameter: To 
be able to talk about a path in parameter-space one must assume that Q is a classical variable so that it makes sense 
to demand differentiability of the paths. Strictly speaking, this assumption is violated when we take into account the 
quantum nature of Q itself, because in quantum theory particles do not follow differentiable paths anvmorj^l. and 
taking this quantum nature into account is a highly non-trivial task (see again |54j and especially [56j for elaborations 
on this topic). However, as long as we are interested in a system where the respective energy scales are well separated 
(to -c M) and we are working in a representation where Q acts as a multiplication operator and not as a derivative 
operator the approximation seems to be valid and is confirmed by many experiments in molecular physics. 

Now that we understood the first half of the problem, let us go back to the full problem (|5.3j) and see how we can obtain 
solutions to the eigenvalue problem for the full Hamilton operator H acting on the fullg space 7i: Assuming that we 



16 To get some intuition for this problem one could think of a quantum particle under the influence of an external magnetic field Q. 
Assume that we can solve the Schrocdinger equation for each fixed value of the magnetic field. Then the question would be: How does the 
quantum system describing the particle change when we manually change the external magnetic field? 

17 The argument turns out to work out in more general situations, see for example 1541 . 

18 In fact, the path integral measure is concentrated on paths that are arbitrarily discontinuous. 
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have solved the eigenvalue problem for H{Q) with Q as an external parameter and we already know its cigenfunctions 
Xi(<?j Q) we start with the Ansatz 



Plugging this into (|5.3[) leads to 



H ^ = J2 (- Mi + A* (Q))& (Q)Xi (?5 Q) = A Q E & (Q)Xi (Q\ 1) 



(5.6) 



(5.7) 



Now we multiply by Xk from the left and take the scalar product in Ti. q to receive an equation for the Q-depcndcnt 
coefficients £f 

E ((xu I (-^ A q) I X*) + W)***) C = A Q 4 Q (5.8) 
where both % and £ are functions of Q and the Laplacian Aq acts on everything on its right. This can be rewritten as 

E (E ^(- ihS ^ 9 Q - A kl )(-iM u d Q - An) + Xi(Q)5 ik j £? = A a ^ (5.9) 



when using the connection 



A k i ■= ih{xk I d Q I xi). 



(5.10) 



Thus we see that the effect of the light variables on the quantum-dynamics of the heavy system is twofold: First, 
there is an effective potential Xi(Q)Sik which can be interpreted as the heavy variables experiencing the presence of 
the light variables only through an average. Second, the light variables effectively curve the space as seen by the 
heavy variables, because the ordinary momentum operator Pki = —ihSkidn gets replaced by the covariant momentum 



operator Pki 



-ih{8. 



hit 



jrAki). Using the covariant momentum operators and defining Pf k := YlPgPik the 



effective Schrocdingcr equation for the Q-dependent coefficients reads 



E 



ik 

2M 



Xi(Q)5 ik 



(5-11) 



So far, apart from the assumption that Q can be treated as a classical external parameter (i.e. showing no quantum 
behaviour as long as the dynamics of the g-variables is concerned), everything was exact. The next step would be to 
solve (|5.11|) and thus to compute the coefficients £f in the ansatz (|5.6p . Thus we would arrive at a complete solution 
^"(QjQ) of the quantum system. 

However, for general potentials V(Q, q) (|5.1ip turns out to be too complicated to be solved exactly, thus one has to rely 
on perturbation theory. The easiest route one can follow, which was originally proposed by Born and Oppcnhcimcr 
in the context of molecule physics [57| and carries their name in the literature is to simply approximate 



Pik « Pfa 

which diagonalises the operator in (|5.11[) and leads to simply 



^ + A fe (Q) 



(5.12) 



(5.13) 



This means that the influence of the g-variables onto the quantum dynamics in the Q-sector is just taken into account 
via the effective potential term Afc(Q). One can show that this approximation is justified whenever |A, — \j\ ^> |A Q — A^| 
(see for example |56j ) and better approximations can be obtained by taking into account off-diagonal matrix elements 
or not setting the connection A^i to zero. 

As we said before for the Born-Oppenheimer approximation to work one must assume that the Q-variables can be 
treated classically when we are only interested in the quantum dynamics of the g-system. This assumption is surely 
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justified when working in a representation where Q acts as a multiplication operator and furthermore the characteristic 
energy scales are well separated (|A^ — Xj\ \A a — A^|): To illustrate such a situation with a well known example 
consider a molecule with a number of nuclei with masses M as Q-variables and a number of electrons with masses m 
as g-variables. Semiclassically speaking, the electrons move much faster than the nuclei, therefore when treating the 
dynamics in the electron-sector one can safely treat the slow nuclei as classical variables. 

There have been some efforts to generalise this picture and take into account the quantum nature of the Q-variables 
when treating the dynamics in the g-sector (see |54j ) but the methods described above do not easily generalise to 
this setting. At least to our knowledge there are no results concerning the applicability of the Born-Oppenheimer 
method for non-commuting Q-variables and a generalisation into that direction seems to require some substantially 
new input. 

The difficulties with the Born-Oppenheimer method for non-commuting variables was one of the main reasons to not 
choose the Ashtekar-Lewandowski-representation for gravity but to work with the one we constructed in this article. 
Here our geometrical variables (given by the cotriads) are represented as ordinary multiplication operators and thus 
the Born-Oppenheimer approximation is directly applicable. 



5.1.2 Application to gravity 

We want to employ an approximation scheme of the type described above to analyse the quantum theory for gravity 
plus matter which we have constructed in section [3] Matter couples to gravity only via the cotriads, not via their 
canonically conjugate momenta, and due to the fact that we have chosen to quantise the cotriads as commuting 
multiplication operators a Born-Oppenheimer approximation scheme can directly be applied. 

We note that the system gravity plus matter has two widely separated energy scales: The matter part of the Hamiltonian 
carries information about the mass of the scalar field, which defines a length scale On the other hand, the 

Planck length £p enters the definition of the gravitational part of the Hamiltonian. For typical values of m the quotient 
« 1. In fact this quotient will be much smaller than the quotient -fjf— in molecule physics, so one should expect 
that the Born-Oppenheimer approximation gives results of very good precision. 

One can interpret this large separation of energy scales in a very intuitive way: During typical interaction processes 
between matter particles the geometry of spacetime changes very slowly. Of course, this will not be the case in the 
deep Planck regime, when considering processes with energies comparable to the Planck energy such that the quantum 
nature of spacetime itself has to be taken into account. But for all particle interactions in the semiclassical regime 
with centre of mass energies much lower than the Planck energy this should be a good approximation. Especially 
when considering a wave function that approximates a classical geometry (like the coherent states for the gravitational 
sector considered in section 2]) the influence of quantum (gravitational) fluctuations will be small. So in a sense this 
method allows us to compute the QFT on CS limit of quantum gravity. However, the method we will describe below 
allows us to go beyond this approximation and at least in principle compute quantum gravitational corrections to 
QFT on CS. 

Let us now explain the Born-Oppenheimer approximation adapted to our gravitational theory: The role of the heavy 
variables is played by the cotriads E\(v) and the role of the light variables is played by the matter fields $(u). To be 
more precise we are considering the Hilbert space TL rs = ® V TL V = ® v (li.® (£> Hf), where H.% = L 2 (IR 9 ), Ttf, = L 2 (R). 
E\(v) acts as a multiplication operator on TC®, -P/(w) := i^pd E i^ as a derivative operator thereon, <&(v) acts as 

a multiplication operator on Ti<$ and Tl(v) = it\B&r v \ as a derivation operator thereon. The Hamiltonian is given 
approximatehf^l given by 



H(P,E,n,$) « [C° in (P > E) + C° ot (E)\ ®id*+C*(E,Il,*) , (5.14) 

where (>(£,n,i>) = c£ in (E,tt) + C* ot (E,$) factorizes into products of the form O 9 ® with 6 9 G C{H%) and 
G CiT-Ct). One difference to the Hamiltonian in the toy model described above is that the kinetic terms CL n 
and C^in are not functions of the momenta alone but are functions of the cotriad-operators E]{v) as well. This does 
not spoil the general argument made above, however, the nice expressions in terms of a gauge potential A/~i and its 
associated covariant derivative is not sufficient anymore and there will be extra terms. 
We want to solve the eigenvalue problem 



H(P,E,n,<5>)y a (E,$) = A a ^ a {E,<5>) , (5.15) 



19 In this approximation we keep only the first term in the power series H := C (l + O I — — Zf Cb I I as explained earlier. 
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and for illustrative purposes we will assume that H has discrete spectrum such that ^f a (E, $) are proper cigenfunctions 
and A Q their respective eigenvalues. At least as far as the matter part of the Hamiltonian is concerned, this will be the 
case whenever the spatial manifold of the continuum theory which arises as an approximation to the discrete theory 
is compact without boundary. The non-compact case has to be treated with more care and one needs to perform a 
direct integral decomposition of H.®- However, these technical details will not be essential for our method and we will 
restrict ourselves to the case of compact spatial manifolds without boundary. 
The strategy is exactly the same as in the toy model discussed before: We start with the Ansatz 



(5.16) 



for states \1/ G H.® and assume that we already have some knowledge about the partial eigenvalue problem 

C'l'($,fl-,E)xi($;E) = \ i (E)xi($;E) . (5.17) 

Here the dependence on E has again to be understood in a parametric sense, i.e. we keep E fixed and regard Xi{®\ E) 
as an element of 7i^. Thus, because C"^($,n; E) is just the standard discrete matter Hamiltonian, this amounts to 
solving a certain type of lattice-QFT for an arbitrary (but fixed) background. Of course, in general this cannot be 
accomplished analytically and one will have to employ further approximation techniques, but the interesting point 
is that QFT on a curved spacetime emerges in the analysis of a theory of quantum gravity in the same way as the 
theory of quantised electrons in a given classical potential emerges out of the full quantum treatment of a molecule. 
For technical reasons one has to assume that Xi(^i E) depends on E in an at least twice differentiable manner. 
In order to get the full wave functions ty a (E, $) we follow the Born-Oppenheimer approach as discussed for the toy 
model above, so we want to compute the coefficients ^f(E). This can be done by starting with (|5.15|) . multiplying by 
Xfc(<E>; E) from the left and taking the scalar product in Tt^. Assuming that Xi(^> E) is an orthonormal basis of 
we get 



Xk($;E) 



cL(p,e) 



C a pot (E) 



Xi(*;E)) £?(E)+X k (E)$(E) 



(5.18) 



If we define CjL 

as 



Xk(^E)\ C° in (P,E) 



C 9 pot (E) 



Xi(&\ E) ) as an operator in C(H®) then we can write f|5 . 18[) 



Cl i + 6 ik \ i (E)W(E)=A a $(E) 



(5.19) 



which is a coupled system of eigenvalue equations in Ti^ and its solutions are the appropriate coefficients in the 
expansion (|5.16p . 

So far everything is exact and one could obtain the full set of solutions by first solving (|5.17[) and then use x% an d 
Xi to compute the coefficients in the expansion (|5. . However, even in our toy model we saw that there is in 
general no chance of obtaining analytic solutions. Quantum gravity is surely more difficult than a single molecule, 
so it would be surprising if one could make progress in the full theory without further simplifications. Thus, the 
crudest approximation is to use again the Born-Oppenheimer approximation which in this case amounts to neglect- 
ing all the off-diagonal terms in the matrix C^, that is we neglect the action of P on Xi- Then we can write 



C 9 pot (E) 



Thus we are left with 



Xk {^E) | C°JP,E) + C 9 pot (E) | XifrE)) , * (Xk(^E) | x^;^ C°.JP, E) + C 9 ot (E) 



pot V 



C° in (P, E) + C g pot (E) 



cL(p,e)- 



(5.20) 



In the toy model we saw that using the pure Born-Oppenheimer approximation means that the influence of the light 
variables onto the quantum dynamics of the heavy variables is only effectively taken into account via their eigenvalues. 
This is also the case here: If we would set X l (E) = then we would describe vacuum gravity and our solutions in the 
gravitational sector would not know anything about the matter content of the theory. One can regard (|5.20p as a first 
step towards quantum gravitational solutions which do not neglect the presence of matter fields, however matter is 
only 'effectively' taken into account and not with respect to full dynamics. In a sense this is the quantum analog of 
the semiclassical Einstein equations = (?),„) ^ where the expectation value of the stress energy tensor is used as 



2G 



a source in Einstein's equations. 

At least in principle the strategy to solve dynamics in the fully quantum gravitational setup of the theory defined above 
would be the following: 

1. As a first step consider the Schrodinger equation (|5.17| for the matter part. This can be interpreted as a funda- 
mentally discrete quantum theory for matter degrees of freedom on a fixed background whose continuum limit 
would be quantum field theory (QFT) on a curved spacetime (CS). However, in the discrete setup we are not 
restricted to smooth metrics since the differential operators have been replaced by difference operators. Hence, 

is a well defined operator for arbitrary values of E\(v) as long as it is non degenerate^. Thus, in a sense 
this can be seen as a generalisation of QFT from spacctimcs with a curved smooth metric to "spacetimes" with 
a curved and even discontinuous metric. Nevertheless, if we consider Ej(v) which approximate a smooth metric 
then (|5.17|) is a 'lattice version' of QFT on CS. Note however, that the interpretation is quite opposite to the one 
usually employed for lattice theories: This theory of quantum gravity is fundamentally discrete and a continuum 
theory should be regarded as an approximation which is good for a certain set of states (the scmiclassical ones) 
and bad for other states. 

2. After analysing the matter dynamics on a fixed background one would obtain eigenfunctions Xii&'jE) an d 
eigenvalues \i(E) which are simply functions of E, i.e. it is not enough to consider QFT on CS for a single 
fixed background, but one would need knowledge about the whole class of quantum field theories on an arbitrary 
background. 

Using these functions we could go ahead and analyse the gravitational sector of the theory, i.e. we would use 
equation (|5. 19|) to compute the coefficients £f in ansatz (|5. 16|) and thus find solutions to the full quantum theory 
of gravity plus matter. However, from a practical point of view, there is no chance that one can exactly solve such 
a complicated system of coupled differential equations as (|5.19[) . so we must resort to approximation methods 
such as the pure Born-Oppenheimer approximation as described above. Thus, one would use equation (|5 .201) 
to calculate the coefficients In this approximation the only imprint that matter leaves on the gravitational 
sector is the term proportional to Ai(E) in equation (|5.20p . 

However, from a practical point of view the situation is less clear: Even if we would succeed to somehow compute 
the eigenfunctions of the matter Hamiltonian Xi(^> E) (either numerically or using perturbative approaches such as 
an expansion in terms of Feynman diagrams) there still remains one big problem: In step [2] of the scheme described 
above we would need to solve the eigenvalue problem for the operator 

Cl in {P,E)+C 9 vot {E) + y{E) (5.21) 

in order to get the correct coefficients This operator is tremendously complicated, as we have seen in (|3.12|) and 
even setting Xi(E) to zero docs not substantially simplify the problem. For Xi(E) = (|5.2ip is equivalent to the 
Hamiltonian constraint of a pure gravity theory on an algebraic lattice, and there is no hope that one can gain enough 
information about the space of its eigenfunctions in an analytic way. This is the reason why it is so complicated 
to incorporate the matter influence on the gravitational sector in the full theory on a more than formal level. Even 
considering \{E) just as a small perturbation does not help much: Without any knowledge about the spectrum of 
the unperturbed operator it is not meaningful to start a perturbative analysis. 

At this point we want to stress that from a conceptual point of view there are no obstacles in the scheme we proposed 
above, it's just that the operator (|5.21[) is too complicated without further simplifications. One could for example 
employ numerical methods to approximately solve (|5.20|) in certain situations. A second option is to further exploit 
the fact that we want to describe situations where geometry can be treated almost classically, thus it makes sense to 
use a semiclassical approximation, which we will describe in the next paragraph. 

5.1.3 Born Oppenheimer in the semiclassical regime 

What we are interested in is not the complete solution to the spectral problem. It would be sufficient to know the 
Xi(E) in the "neighbourhood" of some three geometry Eq and seek for coefficients £i(E) that at the same time solve 
the eigenvalue equation, say in the pure Born - Oppenheimer approximation, and in addition die off sufficiently quickly 
away from Eq. The problem is of course 1. to actually compute the Xi(E) in a suitable neighbourhood of some exactly 
solusolvableable Eq with sufficient analytical control and 2. to find the eigenvectors 

20 In the degenerate case, QFT on CS methods cannot be used. By construction of our domain of definition, this case does not arise 
when restricting the action of the Hamiltonian operator to that domain. 
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As we are unable at present to say much about the solution of either of these problems we must resort here to a 
very crude approximation: 

Pick any solvable E$ (say Minkowski, FRW, ....) and choose as Xi(E) a coherent state ^ z t from the gravitational 
sector alone peaked on a classical phase space point Z 3 j(v) := - ^7= (Ejq ~ i) ( see section Notice that we pick 
^St(Zi) = Eq independent of i but ^s(Zi) = Pi is not specified at the moment. These states are kinematical coherent 
states, not dynamical ones, which means that they are not granted to be stable under the dynamics for a very long 
time. However, when considering particle interactions the typical timescales on which such interactions happen are 
rather short, and at this stage we want to assume that the gravitational coherent states keep sharply peaked on a 
classical trajectory during these time intervals. This assumption seems natural from a particle physics point of view. 
Clearly, the coherent states \P z are not exact eigenstates of the Hamiltonian operator but at least approximately so. 
Moreover, they are orthogonal to a very good approximation whenever the Zi lie in different quantum cells of the 
phase space. As we will see, this condition will hold automatically whenever the Xi(Eq) differ sufficiently from each 
other. We now make the following Ansatz for the "eigenstate" of the physical Hamiltonian 

9(E,<j>)=Y i Ci y Zt (E) Xl (frE ) (5.22) 

i 

for certain complex numbers a (the ^ Zi themselves are normalised). Notice that we have frozen the E dependence of 
Xi{4>] E) at Eq so that C 9 actually does not act on it. Therefore the Ansatz (|5.22p is not really of Born - Oppcnheimcr 
type. Inserting into the eigenvalue equation we obtain the exact equation 

H^ = Y, Cl [(C 9 + X l (E ))^z l KE) Xl (^E )=A^ (5.23) 

i 

and thus due to the orthonormality of the Xi 

d {(C s + X^Eo) - A.)* Zt ](E) = (5.24) 

for all i. Taking the inner product with ^lz 3 we see that, due to our assumption on the Xi(Eo) the equation for j =/= i 
is trivially satisfied (approximately) while for i = j, c, 7^ we obtain 

< Vzv&Vz, >~C 9 (£o,P*) = -Xi(E )+A (5.25) 

which we use as a condition on the Pi which therefore will mutually different from each other as long as the Aj(-Eo) 
are. Notice that our assumption < ipzn i'Zj >~ $ij is therefore self - consistent. 

However, as already stated, this procedure is unsatisfactory for several reasons of which we mention two: First it 
does not really follow the Born - Oppcnheimcr spirit and does not respect the true interaction between geometry and 
matter. This is also obvious from the fact that we nowhere used that the geometry operators were commuting, what 
we did could have been done in LQG as well and has already been done |16j . The only new ingredient is that we 
took into account the backreaction in the sense of expectation values. Second, the above method can never reveal the 
spectral values of the physical Hamiltonian, because we can always choose the Pi to obey (|5.25|) for any choice of A. 

Suffice it to say that a lot of work has to be invested in order to make the Born Oppcnheimcr Ansatz unfold its 
true power. To get more insight into the problem, we consider a truncation of the gravitational phase space (FRW 
spacetimes) in section [5~3l 

5.2 Construction of a Fock space 

Now we want to describe how one can in general construct a Fock space for a linear (i.e. the Hamiltonian is quadratic in 
the dynamical variables) classical theory. As explained above this situation emerges as step[T]in the Born-Oppcnhcimcr 
analysis of our theory: When solving (|5.17[) and only considering matter dynamics the cotriads can be regarded as 
fixed non-dynamical functions because we have chosen a representation where they act as multiplication operators in 
the full theory. We follow largely [TB] and then apply this scheme to the present situation. We assume the reader 
is familiar with the basic definitions concerning Fock spaces, for a more detailed review on the construction of Fock 
spaces in the context of QFT on CS we refer the reader to [55] , 

Let Tii be a Hilbert space (usually called the one particle Hilbert space) with inner product (• | ip = [tt, <j)\ a basis 
of Hi and Q,V £ £(Hi) linear symmetric operators thereon (i.e. Q) = Q,V^ = V). Assume further that their 
inverses Qr x , "P -1 exist on Tii- The Hamiltonian H of a linear classical dynamical system can then be interpreted as 
a functional 

H : Hi -> R; ^ = [tt, c/>] ^ H{^) = 1 (tt | V \ tt) 1 + \ {<f> \ Q | 4>) 1 (5.26) 
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for given operators Q,.,V . 

Further assume that n, 4> form a canonical pair, that is, the Hilbert space Hi has also the structure of a symplectic 
manifold^! (if dim(7ii) is finite and label its basis elements we have {tti, fyi} = Sw and {717,717'} = {</>;, </V} = 
). We are now aiming for a 'quantisation' of that classical linear dynamical system, namely we want to construct a 
Hilbert space H which carries a representation of the classical observable algebra given by the Poisson algebra of 7r 
and <p. 

One convenient way of doing this, which naturally allows for an interpretation in terms of 'particles' is the construction 
of a Fock space: That is, we are implementing a unitary map U on H\ and use this map to construct complex functions 
z(tt, (j)) such that {z;, zi>} = idw where z is the complex conjugate to z. But this is just the commutator algebra of 
creation and annihilation operators and it is well known that such an algebra can naturally be represented on J- s (Hi), 
the symmetric Fock space over Hi. 

Let us start with the Hilbert space Hi and complexify it. Now choose a subspace H^ + C Hi such that there exists a 
unitary, real-linear map U :H i— > Hf + . Using this map we can rewrite the Hamiltonian in the form 

H=(z\u\ z) n c+ (5.27) 

where z G T6[ + and z is the complex conjugate to z in T6[ + and uj G C(Hf + ) as defined below, z is obtained by using 
this unitary map U through^! 

z := -^UiV^-iV^n) (5.28) 
v 2 

V := [7?-l/2( 7 ,l/2 Q7 ,l/2 ) l/2 7 ,-l/2jl/2 : (52Q) 



and uj is given through 



uj := UV^QV^U- 1 (5.30) 



Note that V is a symmetric operator with well defined inverse if Q and V are. 
Then one gets 

(z\uj\z) n c + = ^(UiV^ + tV^Tr) I UV- 1 QV- 1 U- 1 U | (V<f>-iV- 1 n)) n ^ 
- (V^ + iV^w | V- 1 QV- 1 {V<t> - iV- 1 Tt)) l 



- ((V-yV<j> | Q4>\ + - (tt I (V-yV- 1 QV- 1 V- 1 n) 1 



2 

1 / 
2 

+ l - (tt I (P- I )tp- 1 Q0) 1 - i ((O-^t^-^tQt^-l)^ | 7r ) i 



+ -( 7 r|I?- 2 Q0) 1 --MP-2Q0) 1 
- \{<P\Q<^)i + \^\'P^)i 

= H (5.31) 

where in the second equality we used that U is a unitary map, in the fourth equality that Q, V and T> are symmetric 
and in the fifth equality that T> is real and T>~ 2 QT>~ 2 = V '. This last equality can be seen by explicitly writing T> in 
terms of Q,"P and then multiplying both sides first by V^ Yj2 and then by (p 1 /2gpi/2)i/2^ 

From equation (|5.28p one can verify that z, z indeed form the Poisson algebra of latter operators and therefore the 
Hamiltonian (|5.26[) can be defined as a linear operator on J 7 s (Hf + ) in the form (|5.27[) . 

In our case the one particle Hilbert space is given by Hi = l2(V(a)), the space of square-summable functions over the 

set of vertices V(a). The inner product on this space is simply given by (■, -)i : HixHi — > C; /, /' i— > f( v )f'( v )- 

«ev{a) 



21 We will omit the details here, but in the field-theoretic case one has to consider smeared versions of n, <j> with appropriately chosen 
test functions to turn Hi into an symplectic manifold. 

22 It will become obvious why the operator T> must have this rather complicated form in the next paragraph. 



29 



Field coordinates are given by tt(v), 4>{v) and these form a canonical pair with the Poisson bracket {-k(v), 4>(v')} = S vv / . 
If we consider the cotriads on each vertex as prescribed functions then the "classical" Hamiltonian is given by 

H = \ ^ det(E)(v) n2{v) + dct (^)^)^(^)(- A + m ^ v ) > ( 5 - 32 ) 

with the discrete Laplace Beltrami operator A given by 

A = 2dct( ^ )( ^ V /[ det ^)(")^ £eJM ^^^)^^)<^)^(^)V7] (5.33) 

with the forward and backward derivative, Vj and Vj, acting on functions / £ l2(V(a)) as Vj f(v) := f(v + I) — f(v) 
and VJ f(v) := /(w) — /(t) — /) respectively. Thus Q and "P are given as linear symmetric operators on Hi as 

2 = ' (5 ' 34) 

7> = dct(£)(v)(-A + m 2 ) . (5.35) 

So we can construct z, z such that they have the desired Poisson brackets and we can represent them as creation and 
annihilation operators on the Fock space !F s (l2(V(a))). 



5.3 A minisuperspace example: QFT on FRW— spacetimes 

We want to illustrate the Born Oppenhcimcr Decomposition as well as its semiclassical approximation explained 
in section 15.1.31 for a case for which we have sufficient mathematical control over the equations. We will choose a 
"hybrid model", that is, homogeneous, isotropic FRW - minisuperspace model coupled to inhomogcncous matter 
which recently was advocated in the context of Gowdy models [5Hj- That is, in the language of the previous section, 
we pick as a spatial geometry a homogeneous cotriad Eq and pick as a neighbourhood around it all homogeneous 
cotriads E. This neighbourhood confines us to the homogeneous sector of GR while matter is treated as a QFT. In 
the classical theory this leads to consistent equations of motion because the equations of motion for geometry only 
depend on the integrated matter fields. 

We will assume that the spatial topology of the universe which arises in the classical limit of this quantum theory 
is given by £ = T 3 . In the classical continuum theory for these model the cotriad is given by 

K{a):=E8 l a (5.36) 

where E is a homogeneous function of time only, related to the scale factor by a = \E\. Its canonical momentum 

P?(a) := P5? (5.37) 

is also homogeneous and, in terms of r-derivatives takes the form P = —AEE. Assuming that these variables emerge 
from some averaging procedure 



T 3 
T 3 



with Vq the coordinate volume of T 3 they have canonical Poisson brackets 



{P,E} = ^- (5.39) 
6Vq 



In the FRW - minisuperspace model the physical Hamiltonian H p h ys reduces to Hph ys = g^r J r3 day/ (C 9 + C^) 2 ((t) = 
(jg _|_ (j4> because the spatial diffeomorphism constraint identically vanishes. Inserting the ansatz for the elementary 
phase space variables in the general expression (|2 . 14[) the gravitational part of the integrated Hamiltonian is simply 

OTA p2 

8k \E\ 
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Later we will consider spacetimes with non-vanishing cosmological constant A ^ 0, in this case the gravitational part 
of the Hamiltonian constraint is given by 

Passing to the quantum theory this means we have operators E\(v) = E5\ and (v) = P8\ which do not depend on 
the algebraic lattice point v, E = E- acts as a multiplication operator and P = H 2 P dE as a derivative operator. As 
in the Wheeler-de Witt theory the Hilbert space of this minisuperspace model is given by Ti = I^QR), the space of 
square intcgrablc functions over the real line. In the absence of matter fields dynamics in the gravitational sector is 
then generated by an operator C 9 A that is obtained from the classical expression by replacing E — > E, P — > P in some 
suitable operator ordering. The solutions e\(E) to the eigenvalue problem 

C 9 A e x (E) = \e x (E) (5.42) 

can be given in terms of some hypcrgcomctric functions. If we choose a symmetric operator ordering the eigenvalue 
problem reads 

W ° e% -E 2 ^e x (E) + ™&\E\° ex (IS) ^e x (E) + X\E\ 3 e x (E) = (5.43) 



8k dE 2 AV ' 8 1 >dE 



hi. 



and the solutions e x (E) can explicitly be written in terms of Kummcr's functions of the first and second kind (see 
pjD] for their properties). For the special case A = 0, A = the solutions reduce to ordinary Besscl functions and in 
|61| it was shown that in this basis one can construct 'wave packet solutions' whose quantum evolution stays sharply 
peaked on a classical trajectory in phase space. Thus, these quantum solutions are close to classical ones not only at 
one instant of time but also under dynamical evolution. We will follow here a somewhat different route. 

On such a FRW-spacetime one can explicitly carry out the construction of a Fock space using the method explained 
in the last section: If we want this theory to describe quantum field theory on an FRW universe with spatial topology 
£ = T 3 , the number of vertices of a should be finite. Assuming further that the number of vertices in each spatial 
direction is N we have |V(a)| = A^ 3 . Let and TL(v) be the scalar field and its canonical momentum at each vertex 
v of the algebraic graph. Then the classical Hamiltonian for a scalar field on a space with cotriad E is given by 



C*:=\ £ U(v)^U(v) + <P(v)\E\ 

veV(a) ' ' 



E 2 



flat 2 

1 ■ m 



$(v) (5.44) 



with the ordinary flat lattice Laplacian Afl at := S IJ \7j\7j . Using spatial homogeneity and flatness one can introduce 
a discrete Fourier transform from V(a) to K(a), the associated momentum space. K(a) is discrete (because V{a) is 
bounded) and bounded (because V(a) is discrete). If we choose each of the dimensions in T 3 to be of unit coordinate 
length and choose to embed V(a) uniformly in T 3 , then the vertices are labeled by triples v = (v\, 1)2,1*3) with 
vi 7 V2,V3 G {jj, jr, . . . , jj}. Thus, for each function f(v) we can define its Fourier transform 

/k:=£/(w)e~ < *"' (5-45) 

V 

with k — (fci, fe, £3) and fei, fe, k% € {2-7T, Air, . . . , 2Nir} =: K(a). The inverse transformation is then given by 

/(") = ^3 E f* e+ik ' v ■ ( 5 - 46 ) 

keK(a) 

Furthermore we have the identities X) ek ' < - v ~~ v ' = 5 v y and -ws^2e^ k ~ k ~ > ' v = <5fc,fc' where the deltas on the right 

k v 

hand side are both Kronccker deltas. 

One can show that the Fourier transform of the flat Laplacian acting on a function f(v) is given by 

f k e +ik ' v , (5.47) 



Aflat/(w) := -^3 Yl 

keK(a) 



-AN 2 Y sin 2 (-^) 
7=1,2,3 



which in the continuum limit just gives the usual expression lim Aa^fk = —k 2 fk- Here we chose a naive discretisation 
for implementing the derivative on the lattice. Improved discretisations can be found in the context of perfect actions 
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|51| that provide a framework in which discretisation artefacts can be avoided. For this article we will restrict 
our discussion to the standard (naive) discretisation procedure and will present the discussion on improved actions 
elsewhere. 

Thus, we can use the Fourier transform and write the Hamiltonian as 



$ fc i£i 3 



k<£K{a) 



AN 2 



V sin 2 (-^-) 



1=1,2,3 



(5.48) 



In Fourier space the symplectic structure is given by {Hk, ®k'} = N 3 6k,-k' where — k = (2irN — k). Introducing the 
Hilbert space TLi := l2(K(a)) with inner product (/ | g) k := 773 X) /fcff-fc we can write 

keK(a) 



= \ (n I v I n) fc 



1 



:* 1 2 1 *> fc 



(5.49) 



with V := ^3 and Q := \E\ 3 



2^ + m 2 



with 7 (fc) := AN ^ sin (jb) both simply acting by multiplication. 

7=1,2,3 

Here we are in the special situation that [P, Q] = as linear operators on Hi, thus we can explicitly construct the 
algebra of creation and annihilation variables: For commuting V and Q we get 



V 



~Q~ 


1/4 


\E\ 3 / 2 




1/4 











and thus 



Zk ■ = 



V2 



|£nl 3/2 



r 7 2 W 
£ 2 



1/4 



<$>k-i\E\- 3 l 2 



7 2 (fc) 
£ 2 



1-1/4 



(5.50) 



(5.51) 



One can easily check that {zk ,Zk'} = i&k,—k' ■ The standard representation of this algebra is as creation and annihilation 

operators on the Fock space J 7 (Hi). 

Using these variables the Hamiltonian is given by 



(z\uj k \ z) k 



(5.52) 



with ujk = y 1 e' 1 + m2 - Compared to standard QFT on Minkowski space, where the spectrum takes values u>k = 

Vk 2 + rn 2 there are two differences: First, we are considering a discrete theory, thus the momentum squared gets 
replaced by 7 2 (fc), a bounded expression that converges against k 2 in the limit N — ► 00. Second, the Zk, ^k and 
therefore the vacuum |0 > depends explicitly on E. 

If we choose this representation and denote by a k , a,k the corresponding ladder operators in the Fock space can be 
characterised through a cyclic vector ('the vacuum') |0) defined through the condition |0) = Vfc. All other 
elements of T can then be constructed by letting the creation operators act on |0) and the Hamiltonian operator 
is given by 



:= ^2 u k a\ a k 

k 



(5.53) 



For the vacuum state |0) the expectation value of the Hamiltonian is simply (0 | H | 0) = Co and if we denote 
n-particle states for particles with momenta k\ , . . . , k n by 



1 10} 



(5.54) 



we get 



H\n) = £> fei ] |n) =:\ n (E) \n> 

i=l 



(5.55) 



'Here the inner product {. | .) is taken on the full Fock space T as opposed to the one-particle inner product {. | .) k 
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where the role of the index n is played by the labels fei, .., k n . Let us now analyse the coupled quantum system gravity 
(in the minisuperspace approximation) plus matter according to the scheme we proposed in section 15.1.31 To ensure 
that we have nontrivial dynamics for the gravity part already in the absence of matter we choose to work with the 
Hamiltonian C A with non- vanishing cosmological constant. 

In order to carry out the Born Oppenheimer programme, we have to construct states ^(E, <I>) = \ n ) C^>^) 

n 

and to compute the coefficients £(a) in the Born-Oppcnheimcr ansatz. ft is convenient to first pass to a new set of 
canonical variables which are adapted to the specific form C 9 A and then to quantise the system in these new variables 
for which coherent states can be more easily constructed. We introduce the following variables 

E := y^Vo-sgn(E)\E\^ , then {P,E} = K (5.56) 



and all other Poisson brackets arc vanishing. Here we neglected terms involving a derivative of the signum function 
since these terms only contributes when E = 0. The quantity P can be defined as long as E ^ which is given in the 
classical theory. In the quantum theory we will modify the corresponding Schwartz functions as before by means of 
the function s e (E). In these variables C A takes the form 

which is almost of form of a harmonic oscillator except that it has the wrong sign in the kinetic term. We now choose 
a representation on L,2(M.,dE) in which P is implemented as —il 2 P d/dE and E is represented as a multiplication 
operator. The quantum Hamiltonian is thus 




d 2 3 



A 8k dE 2 4k v ; 

This motivates to construct coherent states which simply those of that harmonic oscillator which results when switching 
to a positive sign in front of the kinetic term in 15.571 It is then obvious that these coherent states approximate the 
operators P and E semiclassically well i.e. 



(n o ,E^ Zo )=E = ^Wo- S gn(E )E-\ (*| , P#| ) = P = V3Vb^= (5.59) 




-2 

'3 S 

as long as we are not close to a classical singularity Eq = where the function s e is constant, see also discussion in 
section HI Therefore these states will also give the correct classical value for the expectation value of C A . 

The exact Born Oppenheimer programme now consists in the following: We are looking for the (generalised) eigen- 
values fi and the system of corresponding eigenfunctions n t— ► (E) of the equation 

[C 9 A + A„]£„ = fi Xn (5.60) 

where A„ is a multiplication operator. This can now be done by using quantum mechanical perturbation theory: 
Multiply this equation by n and treat kC a as the "free" Hamiltonian and n\ n as a perturbation potential (notice that 
K is an extremely small parameter. Assuming that the spectral problem for C 9 A has been solved we can then compute 
the /i and the corresponding systems £„. This is beyond the scope of the present paper and we plan to come back to 
this problem in a future publication. From the resulting eigenfunctions ^^(E,^) one then has to select those which 
are peaked on a given spatial geometry Eq and a given extrinsic curvature Po- 

The approximate (scmiclassical) programme consists in choosing coherent states Xn '■= ">Pz n with $t(Z n ) = Eq, ^s(Z n ) = 
P n for some given Eq and to satisfy the eigenvalue equation in the sense of expectation values only which results in 
conditions which are approximately given by 

C 3 A (E ,P n ) + X n (E ) = fx (5.61) 

which in this case can be trivially solved for P n . 

6 Conclusions and Outlook 

In this article we investigated the question of how one can understand (matter) QFT on a fixed background starting 
from first principles, i.e. from a theory where gravitational and matter degrees of freedom are treated as quantum 



33 



variables. In contrast to ordinary QFT on (curved) spacetimes such a theory describes quantum matter on a quantum 
spacetime. We performed an analysis using Born-Oppenheimer methods and it turns out that (matter) QFT on a 
fixed background emerges as an intermediate step when trying to find solutions for the full quantum theory. Within 
the Born-Oppenheimer approximation one assumes that the change in spacetime geometry is very small on typical 
timescales of particle interactions. This is obviously the case in situations where one expects ordinary QFT to be valid 
but will be violated in the Planck scale regime where the notion of a smooth background spacetime is not appropriate 
any longer because quantum fluctuations of the quantum geometry can not be neglected even for very short timescales. 
However, within its regime of applicability, this assumption means that the geometrical operators, which encode the 
gravitational degrees of freedom, can be well approximated by their classical counterparts as long as we are only 
interested in matter dynamics. 

Classically, gravity couples via the spatial metric to standard matter (at least for scalar fields and gauge bosons) 
and not via its canonical momentum. In order to be able to apply the Born-Oppcnhcimcr approximation scheme we 
need to require that the operator analog of the spatial metric can be well approximated by the classical metric in 
some regime. This is surely the case when geometrical quantities are represented as multiplication operators, but is 
not obvious when these quantities are represented as derivative operators because the derivative operators might not 
commute among each other. 

In ordinary LQG we are exactly in the latter situation: The flux operators, which are the quantum versions of densitised 
triads smeared along 2-dimensional surfaces, are represented as derivative operators on the Ashtekar-Lewandowski 
Hilbert space 7Yal and furthermore form a non-commuting operator algebra. Thus it is not immediately obvious in 
which sense they can be used to approximate a classical geometry in the sense of the Born-Oppcnhcimcr approximation. 

This led us to consider a new algebra of operators and a representation thereof which is similar to ordinary LQG in 
the sense that its Hilbert space carries a basis of states which arc defined on certain graphs but deviates from LQG in 
the choice of representation: We worked in a representation where cotriads are represented as ordinary multiplication 
operators and their canonical momenta as derivation operators. By this we obtain a theory of quantum gravity defined 
on an algebraic graph to which we can directly apply the Born - Oppenhcimcr methods because the matter degrees of 
freedom now couple to ordinary multiplication operators. Since the quantum theory defined in this article is an AQG - 
inspired quantisation on a fixed algebraic graph where the notion of a manifold and therefore a classical spacetime has 
a priori no interpretation, we need to analyse its semiclassical sector in order to be able to make contact to ordinary 
QFT. 

Using these methods we saw how QFT on a fixed background emerges out of a full quantum theory and moreover 
were able to give a conceptually framework that in principle allows to calculate backreactions from quantum matter 
onto the quantum dynamics within the gravitational sector. 

However, from a practical point of view there is little chance that one will be able to perform these calculations analyt- 
ically as they require a substantial amount of input concerning the quantum dynamics from the gravitational sector. 
This means that in order to be able to fully understand the problem and to be able to calculate corrections to particle 
physics processes from first principles one will need to get more control over the solution space of the gravitational 
Hamiltonian. 

As in previous contributions to this subject, there exists a 'semiclassical route', i.e. one can consider states that can 
be written as a tensor product of certain appropriately defined Fock states for the matter sector on a fixed background 
and coherent states (which are sharply peaked on a classical geometry around this background) for the gravitational 
sectors. However, while backreaction effects can be taken into account in the sense of expectation values, these states 
are not eigenstates of the quantum Hamiltonian, therefore do not shed light on the true spectral problem and in 
particular do not really follow the Born - Oppenhcimcr scheme. 

Our analysis aims at a complete understanding of particle scattering processes within a full theory of quantum gravity 
for which some knowledge on the spectrum of the physical Hamiltonian and in particular its vacuum sector might 
prove useful. In order to illustrate our method we tested the framework using a FRW minisupcrspacc model for this 
quantisation in section 15.31 There we could see more explicitly how the presence of quantum matter fields can be 
taken into account into the dynamics of the gravitational sector using the Born-Oppcnhcimcr method. 
Of course, in principle there is no obstruction to use this method to analyse how a quantum field theory for matter 
emerges out of a quantum gravitational setting as loop quantum cosmology. Contrary to full LQG the geometrical 
operators are represented as multiplication operators in loop quantum cosmology as well and it will be interesting to 
study the relational between our method and the results obtained in [46] . 

It is challenging to generalise the Born Oppenhcimcr method to representations for which the three geometry is repre- 
sented by non commuting operators such as the non-commuting flux operators of LQG. But this requires a completely 

24 LQC is a based on a representation that mimics that of LQC, it is not of Schrodinger type in contrast to the representation used in 
this paper. 
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new input. 
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A Infinite tensor product Hilbert spaces 

Infinite tensor product (ITP) Hilbert spaces emerge as a generalisation of the tensor product of a finite number of 
Hilbert spaces. They were already studied by von Neumann in the thirties and we want to give a brief summary of 
the material that is relevant for our purposes. For a more detailed discussion see for instance [12] or von Neumann's 
original article [62] . 

Let Tii be Hilbert spaces and i G X for some index set X. In general one can allow for arbitrary cardinality of X but 
for our purposes it will be enough to consider \X\ = N. 

The ITP Hilbert space TC® is defined as the closure of the finite linear span of vectors of the form <£)/:= ®ifi for 
fi G TLi with respect to the inner product 

{®f,® l f)H^~^\{fuf^H i ■ (A.l) 

i 

Here the infinite product 1[\ Z{ of complex numbers Zi = \zi \ exp(z^) is defined via 



n 



n 



cxp U'X>i (A.2) 



if the absolute value Y[ \ z i\ an d the phase ^24>i both converge. In this case we call z :~Yl z i convergent. Otherwise 

i i i 

we set Y[ z i = 0- If the absolute value Y[ \ z <\ converges but not necessarily the phase, we say that Yl z i IS quasi 
convergent. One can show that for z = Yl z i 7^ given any 6 > there exists a finite subset X$ C X such that 

i 

\z — ] | < S. We will only be interested in elements <S>/ that have non-vanishing norm. 

ids 

On a.® one can define different notions of equivalence: Two vectors ®/ and (§>/' are said to be strongly equivalent if 
and only if | — 1| converges. Strong equivalence classes of vectors £3>/ will be denoted by [/]. It follows 

that (® f , ®f>)n 9 = if either [/] ^ [/'] or [/] = [/'] and (/*, f<) Hi = for at least one i. 

If we set (z ■ f)i := Zifi then <E> z -f = (II)®/ f & ^ s to hold if n Zi is not convergent. Provided that there exists 

such z we say that <8> / and <£> /< are weakly equivalent if and only if [z ■ /] = [/] . These weak equivalence classes will 
be denoted by (/). One can show that to functions ®/,®/' are in the same weak equivalence class if and only if 
S I K/ij/i)'Hil — 1| converges. Hence, strong equivalence implies weak equivalence. 

One can show that the closure of the span of all vectors in the same strong equivalence class [/], denoted by W^, 
is separable, consisting of the completion of the finite linear span of vectors of the form ®fi where /' = /j for all 
but finitely many i. The whole ITP Hilbert space is the direct sum (in general an infinite sum) of these separable 
subspaces, = X)^m- also T~tf f ) be the closure of the finite linear span of vectors <g>/' with (/') = (/). Then 
[/] 

the strong equivalence subspaces of are unitarily equivalent with corresponding unitary operators of the form 

U z ®f := ®z-f with Y[ z i quasi convergent. 

iex 

We see that although the ITP Hilbert space H.® is not separable and might look "too large" at first sight, it decomposes 
into a direct sum of separable Hilbert spaces H?^ when taking into account strong equivalence classes of vectors [/] . 
In QFT on a fixed curved background the theory is naturally formulated on a separable Hilbert space, one for each 
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background spacetime. So in a sense these ITP Hilbert spaces emerge as a natural generalisation when going from 
QFT on CS to quantum gravity, which should incorporate QFT on all possible backgrounds simultaneously. 
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